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ABSTRACT 

We present 5 to 36 /xm mid-infrared spectra of 82 young stars in the ~ 2 Myr 
old Chamaeleon I star-forming region, obtained with the Spitzer Infrared Spectrograph 
(IRS). We have classified these objects into various evolutionary classes based on their 
spectral energy distributions and the spectral features seen in the IRS spectra. We 
have analyzed the mid-IR spectra of Class II objects in Chamaeleon I in detail, in 
order to study the vertical and radial structure of the protoplanetary disks surrounding 
these stars. We find evidence for substantial dust settling in most protoplanetary disks 
in Chamaeleon I. We have identified several disks with altered radial structures in 
Chamaeleon I, among them transitional disk candidates which have holes or gaps in their 
disks. Analysis of the silicate emission features in the IRS spectra of Class II objects in 
Cha I shows that the dust grains in these disks have undergone significant processing 
(grain growth and crystallization). However, disks with radial holes/gaps appear to 
have relatively unprocessed grains. We further find the crystalline dust content in the 
inner (< 1-2 AU) and the intermediate (< 10 AU) regions of the protoplanetary disks to 
be tightly correlated. We also investigate the effects of accretion and stellar multiplicity 
on the disk structure and dust properties. Finally, we compare the observed properties 
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of protoplanetary disks in Cha I with those in shghtly younger Taurus and Ophiuchus 
regions and discuss the effects of disk evolution in the first 1-2 Myr. 

Subject headings: circumstellar matter — infrared: stars — planetary systems: proto- 
planetary disks — stars: pre-main-sequence 



1. Introduction 



Planetary systems are formed out of protoplanetary disks surrounding young stars. Under- 
standing the processes responsible for, and the timescales associated with, the dissipation of the 
disks and the formation of planetary systems is an outstanding problem in astronomy. Most solar- 
mass stars younger th an ^ 1 Myr appear to harbor disks around them; by ^ 3 - 5 Myr most of 



them shed their disks (iHaisch et al 



2001 



Hillenbrand 



2005 



Hernandez et al.ll2008l ). Therefore, de- 



tailed studies of the structure and properties of the 1-2 Myr old protoplanetary disks are critically 
important to our understanding of the key processes governing disk evolution and planet formation. 

Protoplanetary disks are natural by-products of the star formation process, which begins when 
a slowly rotating cloud core collapses to form a central protostar surrounded by an accretion disk 
and an overlying, infalling envelope. The materi al from the envelope "rains dow n" to the disk, 
which then gets accreted onto the central star (e.g. IShu et al.lll987l : lHartmannlll998l ). The observed 
spectral energy distribution (SED) of such systeins has a rising conti nuum in the infrared, and 



they have been classified as Class I sources (jLada 



1987 



Wilking)ll989l ). The envelope eventually 



dissipates by draining onto the disk and/or is dispersed by the outflow/wind from the star-disk 
system, leaving behind a pre-main sequence star surrounded by an accretion disk. SEDs of such 
objects show some what fiat or decreasing con tinuum at infrared wavelengths and they are classified 
as Class II objects (jLadalll987l : IWilking]|l989l ) . It is in these disks surrounding Class II objects that 
planetary systems form. 

The planet formation process begins with the subAtm-sized grains iii the disks stickirig tog ether 
to grow to larger mm- and cm-sized particles (e.g. IWeidenschillinglll980l : iBlum Wurmll2008l ). As 
they grow, the larger grains sink down to the disk mid [-plane; sedimentation of the dust can cause 
significant changes in the vertical structure of the disk ([D'Alessio et al.ll2006l ) . Along with the grain 
growth, m ineralization of the init i ally amorphous dus t grains also take place in protoplanetary 



disks (e.g. ICampins Sz Rvan 



198S 



Malfait et al.l Il998l ) . The larger grains that have settled to 



the disk mid-plane further grow i nto km-sized planetes imals, which through collisional growth, 
eventually form protoplanets (e.g IWeidenschillingj l2008l ). Once these 'planetary embryos' have 
become sufficiently massive (~ 10 Mq), they accrete gas from the disks to form giant planets (e.g. 
Pollack et al.lll996l ). A Jupiter-like gas giant formed i n the disk can gravitationally alter t he disk 
structure by forming radial gaps and holes in them (iRice et al.l l2003l : iQuillen et al.l l2004l ) . The 
disks eventually dissipate either by planet formation using up the material and by various other 
disk dispersal processes such as accretion onto the central star, photoevaporation and magneto- 
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rotational instability induced disk clearing (e.g. IClarke et al.l l200ll : [Alexander &: Armitagd 120071 : 
Chiang &: Murrav-Clavl 120071 ) . Pre- main sequence stars which have dissipated most of their disk 
material sh ow margina l or no excess e mission above their stellar photospheres and are called Class 
III objects (|Ladalll987l : [wilkintjllQSgl 'l. 



The effects of disk evolution are best studied at mid-IR wavelengths. The mid-IR emission 
from a Class II object probes the planet formin g region (~ 1-10 AU) of the disk surrounding it 



(e.g. iDullemond et al.l 120071 : iD'Alessio et al.ll2006l). Dust sedimentation makes t he disks flatter, i.e 



reduce the degree of flaring (jPullemond h Dominikll2004l : ID'Alessio et al.ll2006l ). The shape of the 



used to characterize the vertical distribution of the dust in the disks i 


Kenvon &; Hartmann 


1987: 


Dullemond & Dominik 


2004: 


D'Alessio et al. 


200fi: 


Meeus et al. 


2001 




Pur Ian et al. 


2006). 


The 



changes in the radial structure of the disks also affect the shape of the mid-IR continuum: tran- 
sitional disk candidates - objects w i th ra d ial gaps or holes withi n their disks - (IStroni et al.lll989i: 



Skrutskie et al 
2o"o7l : 



1990 



Espaillat et al.ll2oiol ) and outwardly truncated disks (iMcClure et al 



Calvet et al.l l2005l : lEspaillat et al.l l2007bl: iNaiita et al 



200 



3; 



Brown et al 



20081') have been identi- 



fied by inspecting the mid-IR part of the SED. Additionally, the silicate emission features centered 
at ~ 10 and ~ 20 //m generated by the optically thin dust grains in the disk surface layers can 
be used to probe the size and composition of the dust in the disks. The grain growth and crystal- 
lization, i.e. the degre e of dust processing in protoplanetary disks can be quantified usi ng silicate 



emission features (e.g. iBouwman et al 



2001 



van Boekel et al.ll2004l : ISareent et al.ll2009l ^. 



As part of a large Spitzer Infrared Spectrograph Houck et al.l 12004 ) survey of star 

forming regions within ~ 500 pc, our team has obtained mid-IR spectra of several hundred YSOs 



Purlan et al. 


( 


2006^1. 


Purlan et al. 


( 


2008) 



and iMcClure et al.l (|2010l ). Here we present and analyze the IRS spectra of young stars in the 



Chamaeleon I star forming region; a smaller subset of the Class II disks classi fied and iden t ified i n 



this paper along with those in Ophiuchus and Taurus has been analyzed in iPurlan et al.l (|2009l ) 



where a comparative study of the disk properties in these three regions is presented. 

The Chamaeleon I (henceforth Ch a I) cloud, whi ch is one of the nearest star forming regions 
to the Sun, at a distance of 160-165 pc (jLuhmanll2008l ). con tains a relativel y isolated population of 
young stars with modest levels of extinction towards them (jLuhman l2004al). More than 200 young 
stars are known to be members of the Cha I star forming region (jLuhmanI |2008| . and references 
therein). The young stellar population of Cha I has a median ag e of ^ 2 Myr a nd an age spread of 
3-6 Myr and is thus older than that of Taurus and Ophiuchus (iLuhmanI l2008l ) . 



In this paper we present the mid-IR spectra of 82 young stars and brown dwarfs in Cha I, 



^The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded by 
NASA through the Jet Propulsion Laboratory and the Ames Research Center 
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obtained with the Spitzer IRS. Our sample was initially selected from the known members of Cha I 
based on their mid-IR flux measured by the Infrared Astronomical Satellite (IRAS) and Infrared 
Space Ob servatory (ISO) and from the members ide r itified in the deep Ha and near-IR surveys 



of Ch a I (IWhittet et alJ ll987l. Il997l : IPersi et alJl200d : IComeron et alJll999l . l200d : ICambresv et al 

Objects with mid-IR flux F14 3 um or F 



1998; Oasa et al 



19991). 



in the sample. Several low mass stars and brown dwarfs (jLuhman 



9, jim > 10 mJy were included 



2004al ) were also included later 



based on their IRAC and MIPS fl uxes, which were av ailable by then, and employing a flux cutoff 
lira or F24 ^um > 2.5 mJy (see lLuhman et al.l 120081 ). In Section [5] we describe our observations 
and data reduction. In Section [21 we characterize our sample and classify the objects into various 
evolutionary classes based on their spectral energy distributions and the IRS spectra. A detailed 
analysis of the IRS spectra and a comparison of the observed disk properties of the Class II objects 
in Cha I with those predicted by irradiated, accretion disk models are presented in Section [H The 
vertical and radial structure of the disks based on the observed continuum spectral indices are 
discussed in Section 14.11 and the properties of the silicate emission features and dust processing in 
the disks are discussed in Section 14. 2[ In Section [5l we examine the effect of stellar multiplicity on 
disk structure and in Section [6] we compare the disk properties of Classical and Weak-line T Tauri 
stars. Finally, our main results are summarized in Section [71 



2. Observations and data reduction 



We observed 92 young stellar objects, including very low mass stars and brown dwarfs, in 
the Cha I star forming region. All the targets were observed in the IRS staring mode, in which 
a target was successively observed in two nod positions, at ~ 1/3 and 2/3 of the slit length, 
along the spatial direction of the slit. The full mid-IR spectrum from 5 to 38 ^m was obtained by 
using the low spectral resolution (A/AA ~ 60-120) modules of the IRS: Short-L ow (SL; 5.3-14/im) 
and Long-Low (LL; 14-38 //m). For a few bright targets (see I Watson et al.ll2009l ). the high spectral 
resolution (A/AA ~ 600) modules, Short-High (SH; 10-19//m) and Long-High (LH; 19-37 /um), were 
used instead of LL. Each observation immediately follo wed a high-precision pointing peak-up with 
the Pointing Calibration and Reference Sensor (PCRS: iMainzer et al.ll2003l ) yielding a single axis 
pointing uncertainty of 0.28", that is small compared to the pixel sizes (1.8" in SL, 4.8" in LL, 2.3" 
in SH, 5.1" in LH). 

To extract the spectra, we started with the stray light-corrected, dark current-subtracted 
and flat-fielded basic calibrated data (BCD) products from the Spitzer Science Center (SSC) IRS 
data pipeline. The non-flat-fielded, droop products were used for some targets observed using the 
combination of SL, SH, and LH during campaign 29 to enhance the signal-to-noise. The SSC IRS 
pipeline versions S13.2 and S14.0 were used for objects observed in campaigns 20 to 33, and S16.1 
for campaigns 41 and 42 . We first fixed the permanently bad pixels and the so called 'rogue' pixels 
(see lWatson et al.ll2009l ) flagged by the pipeline by interpolation over the neighboring good pixels 
along the dispersion axis (rogue pixels in SH and LH data were fixed for all campaigns; rogue 
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pixels in LL data for campaign 41 and 42 were fixed ; rogue pixels in SL data were not fixed). We 
used the SMART software tool (jHigdon et al.ll2004l ) to extract and calibrate the spectra. For the 
low-resolution data (SL and LL), sky background was removed by subtracting the same orders 
of the spectra obtained at different nod positions before the point source sp ectra were extr acted 
with a variable-wid th column matched to the IRS point-spread function (see Sargent et al.| [2006: 
Watson et al.ll2009l ). For the high resolution data (SH and LH), since the sky emission is negligible, 
a full slit spectral extraction was performed without sky subtraction. The extracted spectra were 
calibrated by multiplying them with the Relative spectral response functions (RSRFs ), constructed 



from the observed spectra of calibration standards and their template spectra (see IWatson et al 



200S 



Sargent et al.ll2009l ). For each object, calibrated spectra at both the nod positions were 



averaged to obtain the final spectrum. Finally, the SH spectra were truncated at 14 /im and all 
the high resolution spectra (SH and LH) were rebinned to the same resolution and sampling as SL 
and LL. 

Special care was required for the spectral extraction of the following objects: Cha Hq!2, Hn lOE, 
ISO 237 and CHXR 30A. The SLl nodi image of Cha Ha2 suffered from an artifact from a bright 
object in the peak-up array, so we used only the nod 2 for the SL spectrum. ISO 237 and CHXR 30A 
had problems with source contamination in the LL slit; however, careful sky subtraction, source 
extraction with a PSF, and scaling between modules generated a reliable version of the 5-36 //m 
spectrum. Hn lOE had scattered light from a nearby bright object entering in the LLl image and 
we could not completely remove this contribution during the extraction, especially at A > 30 fj,m. 
We present the IRS spectrum of Hn lOE here, but do not use it in our analysis. 

Of the 92 sources observed, spectral extraction failed for the followin g 7 objects: HM8 (2M Jl 1025504- 
7721508), ISO 13 (2M Jl 1025579-7724304) (probably a galaxy fsee iLuhman et al.l boosi )). T4 
(SY Cha, 2M J10563044-7711393), 2M J11070768- 7626326, OTS44 (2M J11100934- 7632178), 2M J11114533- 
7636505 and 2M J11122250-7714512. These sources were too faint and could not be located in the 
slit either due to mispointing or due to insufficient exposur e. One of the ob jects, T9, turned out 
be a background late M giant and not a member of Cha I (iLuhmanI l2004al ) . Two of the sources 
observed, 2M J11183572-79355 48 and 2M J11432669- 7804454, are probable members of e Cha as- 
sociation (ILuhman et al.ll2008l ). Thus in this paper, we present IRS spectra of 82 YSOs in Cha I. 
The log of our observations is given in Table [TJ 



3. Characterization of the sample 



The basic data for the Cha I sample are presented in Table[2j We have spectral type information 
for 74 stars out of 82 Cha I members. Most of them (> 90%) are of spectral type KO or later. Of 
the 7 early type stars in the sample, 2 are B type and 5 are G type. Figure [J shows the spectral 
type distribution of stars in the sample which are KO or later. The spectral type distribution of 
our sample peaks at M5/M6 and is quite simila r to the distribution of the more complete sample 
of the known members of Cha I (ILuhmanI l2008l ) : our sample is representative of the young stellar 
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population of the Cha I cloud. A significant fraction of the objects in our sample with IRS spectra 
are very low-mass stars or brown dwarfs: 19 of them have spectral types M5 or later, of which 7 
are likely to be brown dwarfs (M6 or later). 

We classified our sample objects into classical T Tauri stars (CTTS) and weak-line T Tauri 
stars (WTTS) based on their observed equivalent width of Ha emission (W w„) compiled fr o m the 
literature (see Table [2]). We followed the empirical criteria prescribed bv I White Basril (j2003l ) 
which is based on Wna, but also accounts for the spectral type dependence of the 'contrast effect' 
which affects the observed Who- The Wna values used and the CTTS/WTTS (C/W) classification 
of our sample are listed in Table [2j 

In Table [3] we list the observed properties of the companions of the Cha I sources which are 
known to be in multiple systems. We only list the stellar companions whose membership to Cha I 
cloud has been established. For instance, SZ Cha is reported to have two candidate companions, B 
(2M J10581804-771 7197) and C (2M J10581413-7717088) se parated from the primary by 5.2 " and 
12.5 " respectively (jChez et al.l 119971 : iLafreniere et al.l 120081') . However , these candidates have not 
been confirmed as members of Cha I association (jLuhmanl 120071 . l2008l ). Similarly, T51 is reported 
to have a third c omponent T51C (2M J11122740- 7637017) separated by 12.4 " from the pri mary 
(jChez et al.lll997l ). which is not a confirmed member of Cha I association ()Luhmanll2007l . l2008l ) . We 
do not include them in the table. Most of the multiple systems in our sample have companions at 
separations < 3". As the narrowest IRS slit width is 3.6" (SL), the IRS spectra presented here for 
these objects are the composite spectra of the unresolved components. We take this into account 
while interpreting the spectra of these objects. However, for most of these close binaries, one of 
the components is brighter and dominates the mid-IR spectrum. 

We compiled the photometric data for our sample from the literature. For all objects we used 
the J, H, and K .^ measurements frorn the Point Source catalog of the The Two Micron All Sky 
Survey (2MASS; ISkrutskie et al.ll2006l : ICutri et al.ll2003l ). The I-band data w ere obtained from the 
Third Release of the DEep Near-Infrared Survey of the southern sky (DENIS: lEpchtein et al.lll999l ) 
except for two for which there were no DENIS observations. When available, we also compiled the 
IRA C 3.6, 4.5. 5.8, a nd 8. /im and MIPS 24 fim fiux measurements for the objects in our sample 
from iLuhman et al.l (|2008l ). 



Two mid-M type stars in TableEl 2M J11183572- 7935548 and 2M Jl 1432669-780 4454, are more 



likely members of th e e Cha moving group as suggested by their proper motions (see lLuhman et al 



2008 



Luhmanll2008l ). e Cha is older (~ 6 Myr) and closer to us (d ^ 114 pc) than Cha I and its 



meinbers exhibit larger proper motions than those of Cha I (ILuhman et al 



2000 



Frink et al.lll998l ). 



2008 



Mamajek et al 
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3.1. Extinction correction 



Figure [2] shows the distribution of extinction values (Aj) (see Table [2]) of our sample objects in 
Cha I. More than half of the stars have Aj > 0.8 mag. The wavelength dependent extinction law 
towards mole cular clouds is known to diverge froi n that for the diffuse Interstella,r med iu m f ISM) 



(Math 



law of 



1990) for Ay > 3 mag (Aj > 0.8 mag) (IMcChirel l2009l : IChapman et al.ll2009l ). iLuhman 



2004a) has shown that the ex tinction towar ds Cha I does not follow the standard ISM extinction 



Cardelli et al.l (|l989l ) or iMathisI (|l99d ). Moreover, using the standard ISM extinction curve 
to deredden the IRS spectra causes the 10 //m silicate feature to be sharply peaked at ~ 9.7 



/im, especially for objects with large extin ction (IFurlan et al.ll2009l ). For these reasons we adopt 
the composite extinction law prescribed by McClur j ( 20091 ) to deredden the observed broad band 
flux and IRS spectra of Cha I objects. This method ha s been successfully used to deredden th e 
IRS spectra of young stars in the nearby young clusters (jFurlan et al.ll2009l : iMcClure et al.ll2010l ). 
We deredden the broad band photometry and IRS spectra as follows: obje cts with Aj < 0.8 
mag {A-K, < 0-3 mag), are derddened using extinction law from iMathid (|l99d ) with Ry = 5 (see 



Fur Ian et al. 



200S 



McClure et al.ll20ld ) . For objects wi th Aj > 0.8 mag {A-Ks > 0-3 mag), we used 



the appropriate extinction curves from iMcClurd ( 20091 ) to deredden the photometry and the IRS 
spectra. 



3.2. Spectral energy distributions and classification 

We constructed the SEDs of the objects in our sample from the photometry compiled from 
literature and the IRS spectra. Objects for which the Aj values are listed in Tabled the SEDs were 
generated from the dereddened fluxes. The stellar photospheric flux for stars with known spectra l 
types were derived from the colors of main sequence stars tabulated bv lKenvon HartmannI (|l995l ). 
normalized at the dereddened J-band flux of each object. Objects for which spectral types and 
extinction values are not available, we present the observed SEDs. The SEDs of Cha I members 
are shown in Figures [3] - [7] and those of e Cha members are shown in Figure [HI 

We first classifie d the Cha I YSOs into vari o us SED classes base d on the spectral index between 



2.2 yum and 25 ^tcy (jLada 



1987 



Wilking 



198S 



Greene et al.l 11994 ). The spectral index, ^2-25 is 



computed using the relation 



log 



(1) 



from the 2MASS K^-band flux (Ai = 2.159 /im) and the IRS spectrum (A2 = 25 nm). YSOs with 
"-2-25 > 0.3 are classified as Class I sources; objects with 0.3 > 712-25 ^ — 0.3 are classified as 
Flat spectrum sources; sourc es with —0.3 > 712-2 5 > —2.0 as Class II objects; and 712-25 ^ —2.0 
as Class III sources (also see lLuhman et al.ll2008l ). The 712-25 values and the SED classes based on 



them are listed in Table [2j Objects for which Aj values are listed, the observed and the dereddened 
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values of ^2-25 are presented in the table; the SED classes of these objects remain unchanged for 
both these values. For objects with unknown extinction, the SED classes listed in Table [2] are based 
on the observed ^2-25- The SED classification yields 6 Class I sources, 5 Flat spectrum sources, 68 
Class II sources and 3 Class III sources. 

Our SEP classif ication of the Cha I objects are in good agreement with the classification by 
Luhman et al.l (|2008l ) based on their IRAC and MIPS photemetry of these objects. Of the 74 sources 
that are common between the two samples, only 7 show disagreement between the classifications. 
Four of these, T 25, T29, ISO 1 2 6 and T33A, do not have MIPS measurements available and the 
classification by iLuhman et al.l (|2008l ) is based on the 2 to 8 fim slope which differs from our 
classification based on the 712-25 index. For the other three objects, ISO 86, CedllO IRS6 and CI 



25, ou r classification is consistent with that based on the 2 to 24 /im slopes listed in lLuhman et al 



(j2008l ): however, the classification of these objects based on the slope at wavler igths < 8 um 
differ from that based on n2_25 index, so as to make the classification ambiguous (ILuhman et al 



20081 ). Moreover, the extinction towards these objects are uncertain and the spectral slopes are 
computed from the observed SEDs. These objects are further discussed below, where we explore 
the connection between the SED classes and the evolutionary status of YSOs in Cha I. 

The empirical SED classification based on the IR-slope need not al ways correspond t o the 
physical evolutionary stages pr edicted by the star formation theory (e . g. lEvans et al.ll2009l) . As 



pointed out by several authors (IRobitaille et al 



2007 



Crapsi et al.ll2008l : iMcClure et al.ll20inl ;i. the 



spectral index n2-25 is affected by the viewing geometry of the source and/or the line of sight 
extinction, especially in regions of high extinction. A highly extinguished Class II object could 



easily be misclassified as a Flat spectrum or a Class I source (jRobitaille et al.l 120071 : 1 Crapsi et al 
20081 : iMcClure et al.ll2010l ). In order to disentangle the effect of extinction, and to determine if the 



mid-IR emission from Class I and Flat spectrum sources is envelope dominated, we r nake use of 
the "extinction- free" index between 5 and 12 fim (n5_i2) introduced bv IMcClure et al.l (j2010l ). The 
slope of the SEP measured between these wavelengths is shown to remain unaffected by extinction 
(|McClurell2009l : iMcClure et ^boiol ). 



The n5_i2 values for our sample are listed in Table[2j Figure [9] shows the plot between observed 
712-25 index used for SED classification and the "extinction-free" index n^-u- The regions occupied 
by vario us SEP classes are sh own in the Figure. From an analysis of the well studied protostars in 
Taurus. IMcClure et al.l (j2010l ) have shown that for objects with 725-12 > — 0.2, envelope emission 
dominates the mid-IR wavelengths and for objects 715-12 < — 0.2, the mid-IR emission is disk 
dominated. For < — 2.25, the mid-IR emission is mostly photospheric. Below we briefly 

describe the objects in various SED classes and discuss if they are envelope or disk sources. 

Class I Objects: The SEDs of the 6 Class I objects in our sample are presented in Figure [3] 
and their IRS spectra are presented in Figure [TOl Only 4 of them show 715-12 > — 0.2 indicating 
the presence of envelope around them. They are CedllO IRS4, Cha IRN, ISO 192 and T14a. 
Except for T14a, the other three show strong ice absorption features in the IRS spectra: at 6.0 
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fim due to H2O ice, at 6.8 fim due to CH3OH, NH^, and strongly polar H2O ice, and at 15.2 fim 
due to CO2 ice. The 10 //m silicate feature is also strongly in absorption indicating the presence of 
envelope material around them. 

T14a, on the other hand, show featurless mid-IR spectrum. The extinction towards it is 
relatively low {Aj = 0.4 mag) compared to that towards many of the Class II sources in our 
sample. The optical spectru ni of this source sugges ts that it is a classical T Tauri star (see Table [2]) 
with strong forbidden lines (ISchwartz et al .11 19841 ). T14a is unlikely to have an envelope around 
it. It is a 2.35" binary (jHaisch et al.ll2004l ) and the IRS spectrum presented here is a composite of 
both the components. The contribution from the companion is probably responsible for the rising 
SED and its missclassification Class I source. 

Two Class I sources, CedllO IRS6 (ISO 92) and ISO 86 (CEDllO IRSll), have n^^u < -0.2 
indicating that their mid -IR emission is disk d ominated. Both these sources have been classified as 
Flat spectrum sources bv lLuhman et al.l (|2008l ) based on their 3.6-8 ^m and 3.6-24 ^um SEP slopes 



CedllO IRS6 is 1.95" b inary, but the primary appear to dominate the mid-IR flux (|Persi et al 



2OOII : iHaisch et al.ll2006l l. Both CedllO IRS6 and ISO 86 have been reported to have HCO-h (4-3) 
emission, which traces the dense envelope gas; however, the peak intensity of the HCO-I- (4-3) 
emission in both these sourc es are significantly lower than that observed for the envelope sources 



CedllO IRS4 and Cha IRN (Ivan Kempen et al 



200s 



Hiramatsu et al.l l2007|) . CedllO IRS6 is not 



detected at 1.3 mm or at c r n wavelengths, indicati ng weak envelope emission (jHenning et al 



1993 



Lehtinen et al 



2001 



2OO3I : iHiramatsu et al.l 120071 ). Likewise, ISO 86 is not detect ed longward of 



100 um sug gesting that the source does not have dense envelope material around it (ILehtinen et al 



2001 



20031 ). However, both these sources show ice and silicate absorption features (see Figure [TO 



indicating the presence of cold material along the line of sight. They are likely to be disk dom- 
inated sources with low density envelope around them, which are either viewed edge-on and/or 
extinguished by foreground cloud material. 

Flat spectrum sources: The SEDs of the Flat spectrum sources are presented in Figure [4] 
and their IRS spectra are presented in Figure [TTJ All the Flat spectrum sources in Figure [9] fall close 
to the envelope - disk boundary of = — 0.2, indicating that their envelope emission, if any, 

is weak. Their IRS spectra show weak or no ice absorption features; the silicate features in most 
of them are weakly in emission, indicating very little absorbing metrial between the sources and 
the observer. They are likely disks viewed at higher inclination angles and/or disks with remnant 
enevelopes. 

Class II objects: All Class II objects show presence of broad silicate emission features 
at 10 /xm and at 20 ^m. in their IRS spectra. The SEDs of Class II objects in our sample are 
presented in Figure [5] and their IRS spectra are presented in Figures [12] - [T71 The spectra are 
shown in a sequence of decreasing strength of the 10 //m silicate emission feature and decreasing 
mid-IR continuum slope. T33A does not have 712-25 listed in Table [2] as the band flux is an 
upperlimit. But the IRS spectrum shows strong silicate emission at 10 and 20 /xm (Figure [T2]) . 
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and the n5_i2 index is consistent with disk emission. We classify T33A as a Class II s ource. The 



two B-type stars in our samp le, CU Cha and HD 97300 are well known Herbig Be st ars (iThe et al 



1994 



Hillenbrand et al.lll992l ). We classify both of them as Class II objects (also see lLuhman et al 



20081 ) . We present the SEDs and IRS spectra of these Herbig Be stars, but do not include them in 
our analysis. 

Transitional disks: A few of the Class II disks in Cha I have been identified as transitional 



disk candidates bv lKim et al.l (|2009l ). The SEDs of transitional disks are characterised by a deficit of 
flux at wavelengths < 8 //m compared to that of the Class II median, and signific ant excess einission 
comparable to or hig he r than that of C lass II median longward of ~ 13 fim (jCalvet et al.ll2005l : 
Espaillat et"al]l2007al lbl: Ixim et ahlbood). Th ey are protoplanetary disks with dus t depleted inner 



holes or gaps within theni (jStrom et al. 



198S 



Calvet et al. 



2005 



Najita et al. 



2007 



Espaillat et al 



2007bl : iBrown et al.l 120071 ) . In Section 14.1.31 we present a more quantitative criterion based on 
the continuum indices to identify transitional disks from a la. r ge sam ple of Class II disks. We 
confirm the transitional disk candidates reported by iKim et al.l (|2009l ). and identify an additional 
transitional disk candidate, 2M J11241186- 7630425. We reclassify one of the transitional disk 



candidate reported by iKim et al.l (j2009l ). T21, as a Class III object (see below). The SEDs of 
transitional disks in Cha I are shown in Figure [6] and their IRS spectra are shown in Figure [ 



Class III objects: The SEDs of the 3 Class III objects are shown in Figure [71 They have 
n5_i2 < — 2.25 (see Figure [9]) consistent with their emission being photospheric. 

We assesed the completeness of the objects belonging to various SED classes in our sample 
by comparing them with the Cha I YSOs observed and classified with Spitzer IRAC and MIPS 
imaging of the Cha I cloud (see iLuhman et al.ll2008l ). Our sample contains all the Class I and Flat 
sp ectrum sources in Luhman et al.l ( 2008 ) except three, two of which are later than M3. Compared 
to lLuhman et al.l ([2003), our sample of Class II objects is complete for spectral types Ml or earlier 



and 84% complete for spectral types M5 or earlier. Since our sarnple is flux limited, it has only 3 
Class III sources, compared to the 95 identified by [Luhman et al.l (|2008l ). 



e Cha members: Among the e Cha members, 2M J11183572-7935548 appears to be a 
transitional disk candidate and 2M J11432669-7804454 is a Class II object. The SEDs of these 
objects are shown in Figure [H 



Disk structure and evolution 



The mid-IR excess emission observed for the Class II objects arises from the protoplanetary 
disks surrounding them. The shape of the continuum and the dust emission features in the IRS 
spectrum can be used to study the evolution of disk structure and dust properties. In the following 
we carry out a detailed analysis of the SEDs and the IRS spectra of the Class II objects in Cha I 
and compare their observed properties with those predicted by the irradiated, accretion disk models 



(D'Alessio et al. 



2006 



EspaillatI [20091 ). 
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4.1. Continuum spectral indices and disk structure 



In order to study the structure and evolution of Class II disks in Cha I, we quantify the shape 
of their mid-IR continuum by computing various spectral indices from the dereddened 2MASS K^- 
band flux density and IRS spectra. The continuum spectral indices, nx^^Xi^ were computed using 
equation ([T]) between t he wavelengths 2.159 and 5.7 /xm (n^_fi), 5.7 and 13.4 (nfi_i.s) Mm, and 13.4 
and 31.1 //m (ni3_3i) dPurlan et al.lbond . Mm : Iwatson etliDbond : iMcChire et aLlbniol v These 
wavelengths were chosen because they are relatively free of gas and dust emission features. The 
spectral indices computed for Class II, transitional and Class III objects in our sample are listed 
in Table H 

In Figure [T9l we present a plot of the continuum spectral indices ni3_3i and n6-i3 for the Class 
II and Class III objects in Cha I. Transitional disk candidates identified in Section 3.2 are labeled 
separately. Figure [19] shows that the spectral indices of Class III objects are very close to that of 
the stellar photosphere. Most transitional disk candidates separate out from other Class II objects 
in the spectral indices plot. Two Class II disks, T51 and Hn 5, show ni3_3i < -1.33, bluer than that 
for an optically thick, flat disk, suggesting that the disks around them are truncated from outside 
(see lFurlan et al.ll2009l ). Gravitational interaction with a close compa nion (< 1" separation at the 



1993 



Chez et al 



1994 



distance of Cha I) is known to cause outward truncation of the disks (lArtvrnowicz fc Lubow 
McClure et al.ll2008l). T51 is a 1.9^^ binary (ph ysical separation ~ 314 AU) (jReipurth &: Zinnecker 



19971 : iLafreniere et al.l 120081 ) , but it is unlikely that such a wide companion can 



cause disk truncation. Hn 5 is not known to have close companions. 



4-. 1.1. Dust settling in Cha I disks 



Next we compare the observed continuum indices of the Class II disks in Cha I with those 
pr edicted by irrad i ated, accretion disk models. The models used here are based on the formalism 
of iD'Alessio et al.l (j2006l ). where the vertical disk structure is derived self-consis tently. The dust 



in the disk is assumed t o be composed of olivine silicate and graphite grains (jPorschner et al 
I995I: 



Draine Leelll984l ) which are assumed to follow a power law size distribution of the form 
n(a) oc a~^'^. The disk surface layers are populated with mostly small grains 
amaa;=0.25 //m) and the disk midplane with large grains (amm=0.005 /im, 



a. 



mvn~ 



^max ~ 



=0.005 ^m, 
1 mm). The 

models incorporate disk evolutionary effects such as dust settling by parametrizing the depletion of 
the small grains from the surface layers of the disk. The dust depletion factor, e, is defined as the 
ratio of the dust-to-gas mass ratio in the upper layers to that of the interstellar medium. A value 
of e = 1 implies no dust settling and the degree of settling increases with decreasing value of e. For 
e = 0.001, the small grains in the upper layers of the disk are depleted by a factor of 1000 relative 
to the interstellar dust-to-gas mass ratio, with a co r responding inc rease of large grains close to the 
disk midplane (for details see ID'Alessio et al.ll2006l : lEspaillatll2009l ). 



We computed the continuum spectral indices for a grid of disk models for central stellar masses 
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of 0.5 and 0.2 M .p) and for different accretion rates, inclinati on angles and the s ettli ng parameters e 



( Espaillatll2009l ) . These models are the same as those used in lFurlan et al.l (j2009l ) and lMcClure et al, 



(j20ld ). The st ellar masses ch osen for the models are representative of the masses of objects in our 
Cha I sample (jLuhmanl 120071 ). The continuum spectral indices, ng-is and ni3_3i, were computed 
from the model SEDs following the same procedure as for the data. These indices, computed from 
the models, are compared to those evaluated from the IRS spectra of Cha I objects in Figure [20l 
The observed values of ng-is and ni3_3i indices of Class II objects in Cha I are well explained 
by the disk models with settling parameter e of 0.01 and 0.001 for a range of accretion rates and 
inclination angles, which indicates that the surface layers of the disks in these objects are depleted 
of small grains by a factor of 100 to 1000. The disks surrounding Class II objects in Cha I seem to 
have undergone substantial dus t settling. Similar results have also been found for protopla netary 



disks in Taurus and Ophiuchus (jFurlan et al. 



2006; 



Watson et al. 



200S; 



McClure et al.ll20ld ). Even 



though the disks around Class II objects show evidence for significant dust settling, almost all of 
them have values of n6-i3 and ni3_3i indices > —1.33, appropriate for an optically thick, flat 
disk, suggesting that even the most settled disks observed are flared to some degree. Figure [20] 
shows that the disk models, clearly, cannot explain the mid-infrared continuum indices observed 
for transitional disks. They have altered disk structures and the mid- infrared continuum emission 
from them originates in different regions of the disks. 



4- 1-2. Comparison with disks in Taurus 



Substantial dust settling has a l so been observed for the protoplanetary disks in the Taurus star 
forming region (jFurlan et al.l 120061 : IWatson et al.l |2009|) . In Figure [2H we compare the continuum 
spectral indices of low-mass (spectral type GO and later) Class II objects in the Cha I and the Taurus 
star forming regions. The conti nuum indices o f the Class II sources in Taurus were computed from 
their IRS spectra presented in iFurlan et al.l ()2006l ). The transitional disks in both regions are 
not included in the flgure as they have very different disk structures and the interpretation of their 
continuum indices is not the same as that for full disks (see Section STLS]). The distribution of ni3_3i 
(top panel) and n6_i3 (middle panel) indices are very similar for both regions. Two-sided K-S test 
gives the following values for the maximum deviation between the two cumulative distributions D 
and the signiflcance level of the K-S statistic P%: D = 0.09 and P = 94% for ni3_3i and D = 0.1 
and P = 83% for ne_i3, which indicate that the distribution of these two indices are statistically 
the same for objects in Cha I and Taurus. The vertical dust distribution in the disks in these two 
regions a ppear to be similar. Thus protoplanetary disk s in Taurus region, which has a median age of 
~ 1 Myr (iKenyon &: Hartmannlll995l : lHartmannll200ll ). show simil ar degree of du st settling as that 
of those in Cha I, which is slightly older, at a median age of 2 Myr (|Luhmanll2008l ) . Recent analysis 
of the IRS spectra of < 1 Myr Ophiuchus star forming region has shown that the distribution 
of mid-IR continuum indic es of Class II objects in Ophiuchus is statistically not different from 
that for objects in Taurus (jMcClure et al.ll2010l ). This indicates that most protoplanetary disks 
have undergone signiflcant dust settling by < 1 Myr. This is not surprising as theoretical models 
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predic t timescales as short as <C 1 Myr for grain growth a nd dust sedimentation in planet forming 
disks ([Dunemond &: DominikI l2005l : iDominik et al.l 120071 ) . However, protoplanetary disks in the 
older Cha I region show no evidence for higher degree of dust settling compared to the disks in 
younger Taurus and Ophiuchus regions, indicating that the dust settling is not just a monotonic 
function of age. Even within a region, protoplanetary disks show a large range in the observed 
degree of dust settling. This suggests that other processes such as fragmentati on of large grains 
in the midplane and turbulent mixing must be at work (e.g. lOliveira et al.ll2010l ) along with grain 
growth and sedimentation so as to produce the observed statistical distribution of dust settling for 
< 1 — 2 Myr old protoplanetary disks . 

The distribution of n2-6 index in Figure [2T] (bottom panel), however, is quite different for the 
disks in Cha I and Taurus. A K-S test gives D = 0.24 and P = 5% indicating that the distribution 
of n2-6 of Class II (full) disks in Cha I is statistically different from that of Taurus. A significant 
fraction of Taurus objects have redder (more positive) n2_6 values than Cha I objects. To illustrate 
this further, in Figu re [22] we plot the ni.s si index against n2_6 index for Class II objects in Taurus 
and Cha I (also see iMcClure et al.ll2010l ). On average, the observed values of n2_6 index of Cha 
I objects are bluer (more negative) than that for objects in Taurus. The median n2_6 for Cha I 
(n2_6„^^ = -1.61) is lower than the Taurus median (n2_6„^^ = -1.41). This is partly because of 
the difference in the underlying spectral type distribution of the Taurus and the Cha I sample: 
there are many more objects in the M3-M8 spectral type range in the Cha I sample than in the 
Taurus sample resulting in the spectral type distribution peaking around M5 for Cha I and MO 
for Taurus. For Cha I the median n2-6 index for M3-M8 objects (n2_6„g^ = -1.69) is significantly 
bluer than that for the K5-M2 objects {n2-Q^^^ = -1.48) (also see Section [4.3p and this shifts the 
overall median to the bluer side. 



Figure [22] also shows that Cha I lacks objects with extreme red n2-e: as many as 7 Class II 
objects in Taurus have n2_6 ^ -0.5; no objects in Cha I show n2-6 ^ -0.5. Most of the excess 
emission in the 2-6 iim range in Class II objects comes from the innermost parts of the disks - 
from the inner rim at the dust sublimation radius heated directly by the stellar radiation. When 
the accretion rate is high, the added contributi on from the accretion shock on the stellar surface 
increases the net emission from the inner rim (jMuzerolle et al.l 120031 ) which tends to make n2_6 
index flatter (more positive) and when the accretion luminosity is comparable to or higher than 
the intrinsic stellar luminosity, the n2-6 ~ (jD'Alessio et al.ll2009l ). As can be seen from Figure [22] 
(bottom panel), the upper most octile of the n2_6 distribution (12.5% of the points lie above this 
value) for the Taurus sample is significantly redder than that for Cha I. Some of the objects in 
Taurus which have n2_6 ^ ii2-6„pper oct) '^i^-i DP Tau, DR Tau and DC Tau are the highest accretors 
in that region. High accret ors with large n2-fi va lues have also been found among the Class II objects 
in the Ophiuchus region (jMcClure et al.ll2010l ). However, there is clearly a paucity of such high 
accretors in Cha I as indicated by the absence of objects with large n2_6 indices. This is unlikely 
to be a selection effect. The high accretors in Taurus with n2_6 ^ -0.5 have spectral types earlier 
than Ml. As discussed in Section 3.2, the Class II disks in Cha I sample is complete for spectral 
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types Ml or earl ier. Moreover, T Tau ri stars with high accretion rates tend to have higher mid-IR 
excess emission (jD'Alessio et al.ll2006l ). Since our sample is mid-IR flux limited, it is unlikely that 
we have preferentially missed out high accretors in Cha I with spectral types Ml or earlier. The 
lack of high accretors in Cha I in this spectral range is possibly a result of temporal evolution such 
that Cha I which is older, on average, than both Ophiuchus and Taurus appears to have fewer high 
accretors. T his is consistent with t he general finding that accretion rate on average decreases with 
system age (jHartmann et al.lll998l ). 



4-1.3. Transitional and pre-transitional disks 



The structural diversity among the transitional disks has b een demonstrated from the detailed 



modeling of these objects, particularly in the Taurus region (jD'Alessio et al.l l2005l : ICalvet et al 



2005 



Espaillat et al.l l2007bl . l2008l . l2010|). Objects like CoKu Tau4, DM Tau and CM Aur show 



no or at most small continuum excess emission at wavelen gths < 8 /im indicating that their inner 



optically thick disk s have been almost entirely dissipated (jD'Alessio et al.l l2005l : iNagel et al 



201C 



Calvet et al.ll2005l ). The term 'transitional disks' (hereafter TDs) has been specifically used to 
designate such disks with inner holes. Some of these disks may contain small amounts of optically 
thin dust within their inner holes as in the case of GM Aur. The n2-6 index of TDs typically are 
smaller than the lowest octile ('T'2-6io„er oct«e) ^e seen from Figure [22] (top panel). 

In contrast, objects like LkCa 15 and UX TauA, which show a flux deficit at wavelengths 
< 8 jjLvn relative to the the Class II median, have near-IR (< 4 /im) excess emission consistent 
with that from an optically thick disk. These objects have an inner optically thick disk which is 
separated from the outer optically thick disks by an opticall y thin gap. Such disks with radia l 
gaps in them have been called 'pre-transitional' disks (PTDs) (jEspaillat et al.ll2007bl . l2008l . l2010l ). 
with the name implying an evolutionary stage earlier to that of TDs. PTDs show higher excess 
at shorter wavelengths than the TDs and their n2_6 values fall between the lowest octile and the 
median (n2^6we. oc«ie ^ ™2-6pti5 ^ ^2-6^,d,an) (see Figure [22|). Both TDs and PTDs in Taurus 
show ni3_3i > the upper octile (?^l3-3l„pp,, ^^^ie) (jMcClure et ahlboid ). 

Based on the locations of well modeled TDs and PTDs in Taurus in the ni3_3i - n2_6 plot 
(Figure [221 bottom panel), we classify CS Cha, T25, T54, CHXR 22E and 2M J11241186- 7630425 
in Cha I as TDs with inner holes (also see [Kim et al.[ [20091 : [Furlan et al.l [20091 ). All these objects 
have n2_6 < ™2-6,,„,, and ni3_3i > ni^^^^^^^r ocuu- The disk around CS Cha has already 
been modeled and is found to have an inner hole of ~ 43 AU in radius (jEspaillat et al.[[2007a[ ). The 
loc ations of SZ Cha , T35 and T56 in F igure 1221 (bottom panel), would qualify them as PTDs(also 



see 



Kim et al.[ [20091 : [Furlan et al.[ [20091 ) . Their n2-6 values fall between the median and the lower 



octile and their ni3_3i values are above the upper octile. Modeling of the individual SEDs and 
high resolution submm/mm imaging of the TDs and PTDs are required to constrain the detailed 
structure of their disks. 
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4.2. Analysis of the silicate emission features 

4-2.1. Silicate emission and disk structure 

In this section we explore the connection between the disk structure as revealed by the mid-IR 
continuum spectral indices and the s trength of the silicate emission fea.tures i n the observed spectr a 
of Class II disks in Cha I (see also IWatson et alj [20091 : IPurlan et aPbood : iMcClure et aU boiol ). 
To quantify the strength of the silicate emission we computed the equivalent width of the dust 
emission features centered at 10 //m and 20 //m. The equivalent width is defined as 



W 



A2 



,cont 



Ai 



F: 



d\ 



X,cont 



(2) 



and we integrate between Ai = 8 fim and A2 = 13 ^m to obtain the equivalent width of the 10 fim 
feature,Wio and between Ai = 16 fim and A2 = 28 //m to obtain the equivalent width of the 20 
^m feature, W20- We also computed the integrated flux of the 10 and 20 fim features as 



rA2 

F = (Fa - F; 
^Ai 



\,cont. 



dX 



(3) 



and integrating between the same limits as above to obtain Fio and F2o- In order to compute the 
equivalent widths and the integrated flux of the silicate emission features, we flrst flt a continuum to 
the observed spectra as a polynomial of degree 3, 4 or 5 depending on the shape of the continuum. 
The points used for the continuum flt were chosen so that the dust emission features are small 
within those wavelength r anges: 5.61 - 7.94 um . 13.02 - 13.50 ^m, 14.32 - 14.80 fim , 30.16 - 32.19 
fim and 35.07 - 35.92 fim (jWatson et al.ll2009l ). Objects for which the continuum could not be flt 
by a single polynomial, we used a two polynomial fit, one between wavelength interval 5-14 fLm 
and the other between 14 - 36 fim, making sure that the overall continuum is smoothly varying. 
The values of integrated fiux and equivalent widths of the silicate emission features thus obtained 
are tabulated in Table [5j The dominant source of errors in the quantities estimated from the 
observed silicate features is the uncertainty in the determination of the underlying continuum. In 
particular, the equivalent width of the 20 fim feature, which is broader and generally fiatter than 
the 10 fim feature, is much more sensitive to the choice of the continuum. The errors listed in 
Table [5] assumes an uncertainty of 20% in the estimated continuum underlying the 20 fim feature 
and a 10% uncertainty in the continuum for other quantities. 

Figure [23] shows a plot between equivalent width of the 10 ftm silicate feature, Wio against 
the continuum index ni3_3i for the Class II objects in Cha I. Also shown in the Figure, marked 
as the black polygon, is the region defined by the models in the Wio - ni3_3i plane, for irradiated, 
accretion disk models for a range of values of accretion rates, stellar masses, inclination angles and 
settling parameter e (jEspaillati l2009l : iFurlan et al.ll2009l ). It can be seen from the figure that most 
of the Class II objects fall within the region predicted by the models. Not surprisingly, objects 
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identified as having altered radial disk structures - the outwardly truncated disks and TDs and 
PTDs - fall outside the polygon. Only 4 out of the 8 TDs and PTDs, SZ Cha, CS Cha, T56 and 
2M J11241186- 7630425, show the 10 /nm silicate features in emission and only they are shown here 
(see Figure [THj) . The silicate emission in these disks are generated either by the optically thin dust 
grains with in the inner holes /gaps or at the atmosphere of the inne r 'wall' of the outer optically 



thick disk (Calvet et al 



2005 



Espaillat et al. 



2008 



KimeialJl2Q0g|). TDs and PTDs in general 



show relatively high values of Wm and ni.s -ii, i.e ., they are outliers in Wio and/or ni3_3i as can 



be seen from Figure [Ml (see also lFurlan et al.ll2009l ). Six objects, TW Cha, ISO 91, CR Cha, T33A, 
B43 and CV Cha, (gray solid circles in Figure r23|) show enhanced 10 /im silicate emission. They 
have Wio values greater than the upper octile of the Wio distribution (dashed line in Figure [23]) 



of all the Class II objects. Their Wio values are significantly higher than those predicted by the 
models for typical disks; however, they have ni3_3i values within the range allowed by the models. 

Figure [2H compares the Wio values of Class II objects in Cha I with the integrated flux of the 
10 fim feature Fio. The equivalent width of the silicate feature is a measure of the amount of the 
optically thin dust per unit projected area of the optically thick disk whereas the integrated flux 
of the silicate feature is a measure of the amount of optically thin dust. Both Wio and Fio are 
measures of the strength of the 10 fim feature, but Wio expresses it in the units of the underlying 
continuum. For radially continuous 'full' disks both the equivalent width Wio and the integrated 
flux Fio track each other well as shown in Figure [231 TDs and PTDs show moderate to low values 
of Fio indicating that the amount of optically thin dust contributing to the 10 /im emission is 
relatively small; and yet they show large Wio values. TDs and PTDs have holes/gaps in their 
disks which lowers the continuum emission thereby increasing their Wio- The objects identified as 
having enhanced 10 //m emission (solid gray circles in Figure [24l) have Fio values similar to other 
Class II objects in the sample, but show significantly high Wio values. The large Wio values of 
these objects are likely caused by the reduced continuum emission in these objec ts. Most of the 



conti nuum emission underlying the 10 fxm feature comes from < 1-2 AU of the disks (jD'Alessio et al 
20061 ). One possible explanation is that these disks are opening up gaps in the inner regions (< 



1-2 AU) of their optically thick disks, which will lower their continuum emission. If the gaps are 
sufficiently small (a few AU), one will not be able to distinguish these objects from the ' full' disks 
purely based on the mid -IR continuum; however, they will stand out as outliers in Wio ([Espaillat 



200S 



Furlan et al.ll2009l ). Such objects which are possibly opening gaps within them have also bee n 
identified in the Taurus and Ophiuchus star forming regions as outliers in Wio (IFurlan et al.l[2009l ). 



However, detailed modeling of the SEDs and high resolution imaging are required before the disk 
structure of these objects can be confirmed. 



In Figure [25} we plot the equivalent width of the 20 //m feature, W20 against the continuum 
index ni3_3i. Most TDs and PTDs show large W20 values. The holes/gaps in them are sufficiently 
large to lower the continuum emission underlying both the 10 and 20 /im silicate emission features. 
The objects with Wio ^ ^ 10, upper octUe (gray solid circles), do not have such high W20; their W20 
values, though higher than the median, are similar to those of many other Class II disks in the 
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sample. This is consistent with the supposition that gaps in the these disks are small and are at 
smaller disk radii so that the continuum emission underneath the 20 ^um feature in these disks is 
not significantly lowered so as to enhance W2o- 

Leaving aside the TDs/PTDs and the objects with enhanced Wio, Figure ES] shows that, in 
general, the equivalent width of the 10 //m feature, Wio, increases with increasing value of ni3_3i 
index for 'full' disks. A similar positive correlation is also found between W20 and the ni3_3i index 
as shown in Figure [25j The linear (Pearson) correlation coefficient for these relations is r = 0.3 
and the probability that a randomly drawn sample of the same size could show such a correlation, 
is P < 5%, indicating that the correlations are marginally significant. Such correlations between 
the strengths of the silicate features and the disk geometry are expected for radially continuous 
'full disks' from the disk models: when ni3_3i increases (less dust settling in the disk) degree of 
flaring of the disk increases and the location of the silicate emission moves farther out in the disk 
resulting in a larger emitting area; thus the amount optically thin emission from the warm grains 
in the disk surface layers will increase. However, when the degree of flaring increases, the height 
of the irradiation surface also increases, which increases the amount of heating of the disk interior 
and as a result the continuum level goes up. The net effect appears to be an overall increase in the 
equivalent width of the silicate feature a nd the model predictions agree with such an increase (see 



D'Alessio et al.ll2006l : iFurlan et al.ll2009l ). We note here that the spread in stellar masses, accretion 
rates and inclination angles of the objects in our sample would act to weaken the trend between 
equivalent widths of silicate features and ni3_3i. Despite this, marginally significant correlation 
could be found between the strength of the silicate emission and the underlying disk structure; this 
makes the case for it even stronger. 

Next we look at how the changes in the disk structure affects the strength of the 20 /um silicate 
feature relative to that of the 10 /xm feature. In Figure [26] we plot the ratio of the integrated flux 
of the 20 fim silicate emission to that of 10 /xm emission as a function of the ni3_3i index. We find 
a tight correlation between the F20/F10 ratio and the ni3_3i index for the 'full disks'. The Pearson 
correlation coefficient is r = 0.6 and the corresponding probability is P <C 0.05% indicating that the 
correlation is highly significant. The correlation between the integrated flux ratio F20/F10 and the 
continuum index ni3_3i is much tighter than the correlations that we find between Wio ^-nd W20 
with ni3_3i. This indicates that even though both the 10 and the 20 /im feature strength decreases 
with increasing degree of dust sedimentation in the disk, the strength of the 20 /xm silicate emission 
feature drops faster relative to the 10 fim feature. This is because the 20 /im emission arises from 
larger disk radii compared to 10 /im emitting region. As the disks become less flared because of 
sedimentation, the optically thin emitting area decreases; however, the reduction in the emitting 
area at larger disk radii is more than the reduction in the emitting area closer to the star and 



therefore the the 20 fj,m emission is affected more than the 10 /im emission (see iD'Alessio et al 



2006 



EspaillatJl2009l V 



TDs and PTDs separate out from the 'full' disks in F20/F10 - ni3_3i plane. This is not 
surprising as these disks have different radial structures and both the mid-IR continuum and the 
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silicate emission arise in different parts of the disks as compared to 'full' disks. The objects identified 
as having high Wio (gray solid circles) do not follow the trend seen between F20/F10 and ni3_3i 
for "fuir disks. They show smaller F20/F10 ratio than that for objects with similar ni3_3i values. 
This is suggestive of their disk structure being different from that of radially continuous 'full' disks. 



4-2.2. Silicate emission and dust processing in disks 



We further analyzed the silicate features in the mid-IR spectra of Class II objects in Cha I 
to study the amount of dust processing in protoplanetary disks. Figure [23 shows a plot of the 
flux ratio F11.3/F9.8 versus the equivalent width of the 10 /im silicate feature, Wio- The flux ratio 
F11.3/F9.8 has long been considered as a measure of the degree of dust p r ocessing (grain growth and 



crystallization) in t he inner 1-2 AU of t he disk (jvan Boekel et al. 



Honda et al 



2006 



2005 



Kessler-Silacci et al 



2006 



Olofsson et al.l l2009l ). The larger values of the Fii.3/F9,8 ratio indicates the 



presence of highly processed grains and the smaller values imply relatively unprocessed grains. The 
flux densities Fg g and Fii,3 are computed from the IRS spectra by integrating the flux densities 
within an interval of 0.4 fj,m centered at 11.3 and 9.8 fim respectively. The F11.3/F9.8 values are 
listed in TableO One of the TD/PTD objects, SZ Cha is not shown in the Figure [27l because it has 
strong PAH emission at 11.3 fiui which makes reliable estimate of the F11.3/F9.8 difficult. We also 
show TDs and PTDs in the Taurus region (open star symbols) in Figure [27] for comparison. The 
'full' disks (open circles) show a strong negative correlation between the dust processing indicator 
F11.3/F9.8 and WiQ. The Pearson coefficient for the correlation is r = -0.5 with a corresponding 
probability P < 0.05%, indicating a highly significant correlation. Objects with low Wio values 
show evidence for highly processed dust grains and objects with high Wio values show relatively 
unprocessed grains. 

In the last Section we showed that in 'full' disks Wio is a measure of the amount of the optically 
thin dust responsible for the 10 fim silicate emission, and less settled disks (higher degree of flaring), 
in general, have high Wio. The correlation that we find between Fii.3/F9,8 and Wio also suggests 
that more settled disks show higher degree of dust processing (grain growth and crystallization). 

The TDs and PTDs in both Cha I and Taurus have relatively low values of Fii.3/F9,8, indicating 
that the dust grains contributing to the 10 /xm emission in them are relatively unprocessed. Inter- 
estingly, the objects identified as having enhanced 10 /im emission (solid gray symbols) also show 
low F11.3/F9.8 values, indicating that the optically thin dust in them are relatively unprocessed. 
The high Wio values and the pristine dust cor aposition that these objects show is quite similar to 



that exhibited by PTDs in Taurus and Cha I (jWatson et al.ll2009l : iKim et al.ll2009l : ISargent et al 



2OO9I ) and very different from that shown by 'full' disks. This suggests that as in PTDs, the 10 
fim silicate emission in these disks is likely generated by optically thin dust within the gaps in the 
disks. A more detailed discussion on the pristine dust composition of these disks will be presented 
in a forthcoming paper (Manoj et al. 2010, in preparation). 
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As pointed out earlier, both grain growth and crystalhzation in the inner disk (< 1-2 AU) 
affect the Fn.s/Fg.g ratio. In order to disentangle the effect of crystallization from that of the 
grain growth, we make use of the isolated forsterite feature at 33.6 /im, which is an independent 
measure of the degree of crystallinity a nd, in addition, probes the crystalline co ntent of the cooler 
dust farther out (< 10 AU) in the disk (joiofsson et al-lbood : ISargent et ahlboog ). We compute the 
equivalent width of the 33.6 ^ m feature, W3.S, from equation (l2l) by in tegrating the flux between 



Ai = 32 fim and A2 = 35 fim (jWatson et al.l 120091 : iMcClure et al.ll201Cll ): the W33 values are listed 
in Table [SI Figure [25] shows W33 plotted against F11.3/F9.8 for Class II objects in Cha I. W33 and 
Fii.3/F9,8 appear to be tightly correlated; a linear (Pearson) coefficient of r = 0.4 and an associated 
probability of <C 1%, indicates that the correlation is highly significant. The higher the degree of 
dust processing (grain growth and crystallization) in the inner disk, the higher is the crystalline 
content of the outer disk. F11.3/F9.8 and W33 have also been found to be correlated for Class II 
disks in other star forming regions. In Figure [29] we show W33 plotted against F11.3/F9.8 for low 
mass Class II objects in Taurus. These quanti t ies w ere c omputed from t he IR S spectra of Class 
II objects in Taurus presented in iFurlan et al.l (j2006l ) and IWatson et al.l (|2009l ). Again we find a 
tight correlation between Fii,3/F9,8 and W33 with a (Pearson) correlation coefficient r = 0.5 and an 
associated probability of <C 1%. Similar result has also been found for Class II disks in Ophiuchus 
(|McClure et al.lboiol ). 



Although both grain growth and crystallization in the inner disk can affect the F11.3/F9.8 
ratio, the correlation that we find is more likely due to the crystall ine content in the outer and 
inner regions of the disks tracking each other (see I Watson et al.ll2009l ). Since the crystalline grains 
are generally thought to be formed in the i nner (< 1 AU) d i sks, where the temperature req uired for 
crystallization can be reached (> 1000 K) (Grossman 1972 : Gaij 2001 : Wooden et al. 2005 ). higher 
crystalline content of the outer disk should indicate that the crystalline grains are also present in the 
inner disk. The above results suggest that an increase in the inner disk crystallinity is accompanied 
by a similar increase in the outer disk crystal linity in protoplanetary disks surrounding low mass 
stars (jWatson et al.ll2009l : iMcClure et al.ll2010l ). If the crystalline grains are produced only close to 
the central star, then the correlation that we find between the inner and the outer disk crystallinity 
should suggest efficient radial mixing of the dust in the disk at least up to a radius of < 10 AU. 
If, on the ot her hand, the crystall i ne grains are produ ced insitu in the outer parts (< 10 AU) 



of the disks (Harker &: Desch 



2002 



Wooden et al.ll2005l ). then both the inner and the outer disk 



crystallization mechanisms must operate on similar timescales. 

The changes in the Fii.3/F9,8 ratio have at tinies been interpreted as a l most entirely due to 



grain growth in the inner disk (e.g. iPrzvgodda et al.ll2003l : [Honda et al.ll2006l : [Kessler-Silacci et al 
20061 ). However, the correlation that we find between F11.3/F9.8 and W33 would suggest that 
F11.3/F9.8 is a s good or a better indicator for the cry stalline content of the inner disk as it is for 



grain growth (jWatson et al.l 120091 : 1 Sargent et al.l l2009l ) . Since changes in the grain size affect the 
F11.3/F9.8 ratio, our results also imply that the grai n growth and cry stallization in the inner disk 



must occur concurrently, on similar timescales (e.g. [Apai et al.[ [20051 ). Detailed dust composition 
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analysis of the IRS spectra of T Tauri stars in various star forming regions has sh own that both 
large (> 1 fim) and crystalline grains are present in most protoplanetary disks (e.g. I Sargent et al 



20091 ). 



4.3. Median SED of Class II objects 

We constructed median SEDs for the Class II objects in Cha I from 0.9 fim to 36 ^um using 
/, J, H and Ks photometry and IRS spectra in order to study the average disk properties. Two 
separate medians were computed for objects belonging to the spectral type bins of K5-M2 and 
M3-M8 to compare the emission from protoplanetary disks with different central star properties. 
The median was computed for 25 objects with spectral types between K5 to M2 and for 29 objects 
between M3 to M8. Objects which have only upper limits in /, J, H and Kg measurements were 
excluded from the median. The medians for each of these bins were obtained by first normalizing 
the dereddened flux densities at I, H and Ks and IRS spectrum of each object to the dereddened 
J-band flux of that object, and then computing the median at I, J, H and Kg and the IRS 
wavelengths. We also obtained the quartiles for the medians: the upper and lower quartiles define 
the range above and below the median within which 50% of the flux values lie. The median SEDs 
for the Class II objects in the spectral bins of K5-M2 and M3-M8, along with the quartiles, are 
presented in Figure [301 The gray shaded regions in the figure are bound by the lower and upper 
quartiles of the medians. The photospheres plotted in the figure are that for a K7 star for the K5-M2 
median and an d that for a M5 star for the M3-M8 median, both computed from the photospheric 



colors listed in lKenvon &: HartmannI (|l995l ) and normalized to the J-band flux 



The median SEDs computed for the T Tauri stars (K5-M2) and very low-mass stars and brown 
dwarfs (M3-M8) clearly show mid-IR excess emission and 10 and 20 fim silicate emission features. 
The mid-IR continuum spectral indices, ng-is and ni3_3i, of both the medians are very similar: 
n6_i3 = -0.9 for both K5-M2 and M3-M8 medians and ni3_3i = -0.26 for K5-M2 and -0.28 for 
M3-M8 median. However, the nature of excess emission below 8 /xm is very different for these two 
medians; the excess emission in K5-M2 median begins at ~ 2 fim whereas for the M3-M8 median 
significant excess starts to show up only at 4-5 fim. This is because the ratio of the disk emission 
to the stellar photospheric emission scales as T^, and for cooler stars (M3 or later) i t is difficult to 



see the disk excess above photospheric emission at A < 6 fim (jErcolano et al.ll2009l ). In addition, 
we also find the continuum slope at wavelengths < 8 fim for M3-M8 median is much steeper (bluer) 
than that for the K5-M2 median: K5-M2 median has the continuum index n2_6 = —1.42 while 
M3-M8 median has n2_6 = -1.74. 

It is important to recognize that very low-mass stars and brown dwarfs on average have fainter 
and steeper SEDs at wavelengths < 8 ^m, especially when identifying TD and PTD candidates. 
They are identified by virtue of the flux deficit and steeper slope that they show at wavelengths < 
8 ^m compared to the median Class II SED. A very low-mass object or a brown dwarf with normal 
disk could easily be misidentified as TD or PTD candidate if compared to a T Tauri (K5-M2) 
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median (see lErcolano et al.ll2009l ) 



The e ffect of dust sedimentation in protoplanet a ry disks is best seen in the SED at wavelengths 



> 13 fim (jD'Alessio et al. 



2006 



Furlan et al.ll2006l : IWatson et al.ll2009l ). To demonstrate this, in 



Figure [301 where the median SED of K5 - M2 Class II objects in Cha I is presented {top panel), we 
also show synthetic SEDs generated from the accretion disk models for a 0.5 Mq star with accretion 
rates of 10~^, 10~^ and 10~^ MQyr"^ and the settling parameter e of 0.001, 0.01 and 0.1. It can 
be seen that the median SED longward of 13 yum is best described by disk models with settling 
parameter e = 0.01 and 0.001 for accretion rates of 10^^ and 10"*^ M0yr~^; even for an accretion 
rate of 10~^ MQyr"^, most disks have e < 0.1. The median and the quartiles suggest significant 
dust settling in the Class II disks in Cha I: small grains in the upper layers of more than 75% of the 
disks have been depleted by a factor of at least 10; about 50% of the disks have depletion factor > 
100. 

In Figure [30] (bottom panel), we show the synthetic SEDs generated from the accretion disk 
models for a 0.2 Mq star with accretion rates of 10~^, 10~^ and 10"^'' M0yr~^ and the settling 
parameter e of 0.001, 0.01 and 0.1. The observed median and the quartiles at wavelengths > 13 
^m, where the dust sedimentation effects are seen, are best explained by disk models with settling 
parameter values of e = 0.01 and 0.001 for the above accretion rates indicating significant dust 
depletion in the upper layers of these disks. Both the K5-M2 stars and the low-mass objects with 
spectral types M3-M8 seem to have dust depletion factors of 100-1000 in their upper layers, similar 
to what was found for Taurus (jFurlan et al.l [20061 ). 



However, the match between the model predictions and the median SEDs is not satisfactory 
at shorter wavelengths. The accretion disk models for the 0.5 Mq significantly under predicts the 
observed flux of the K5-M2 median at wavelengths < 8 fini. Also, the upper quartile of the M3-M8 
median at wavelengths < 13 fim cannot be explained by the 0.2 Mq models even for an accretion 
rate 10~* MQyr~^. These discrepancies between the mode l s and the observed medians and the 
possible reasons for them are discussed in detail by [Espaillat[ (l2009l ). Here we only discuss behavior 
of the median at wavelengths > 13 /.im which is relevant for the dust sedimentation effects. 



5. Multiplicity and disk structure 



In this Section we discuss the effect of multiplicity on the structure and dust properties of 
the Class II disks in Cha I. Cl ose companions can gravitat ionally perturb and truncate the disks 
around individual components (ILubow Artvmowica[2000l ). For orbital separations of ~ 200 AU, 
companions truncate the disk from outside (for typical disk sizes of ~ 100 AU); for even smaller 
orbital separatio ns ( < a few tens of A U) the disk could be truncated from inside as in the case 
of CoKu Tau/4 ([Ireland &: Kraud [20081 ). An outwardly truncate d disk is likely to sho w a steeper 
mid-IR continuum slope, i.e. ni3_3i less than that for a flat disk (iMcClure et al.l[2008[ ). Similarly, 
disks truncated from inside will have significantly steeper ng-is continuum index. There have also 
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been suggestions that close companio ns could affect the dust properties in protoplanetary disks 



(Meeus et al. 


2003: 


Sterzik et al. 


2004) 



We compared the observed disk properties of the objects in multiple systems to those in single 
systems in order to study the effect of companions on the disk structure. Figure ET] compares 
the distribution of continuum indices of close (< 1" ~ 165 AU) binaries and wide (> 1") binaries 
among Class II objects with those with no known companions. A two-sided K-S test shows that 
the distributions of the n2-6 and ni3_3i indices for the sub-arcsecond binaries are similar to those 
for the wide binaries and single stars. However, the ng-is distribution for the sub-arcsecond 
binaries appears to be different from that for the wide binaries and single stars at a significance 
level of 97%. The ng-ia distribution for the sub-arcsecond binaries peaks at bluer value than 
that for other objects. We also compared the silicate feature properties, Wio and W20, and dust 
processing indicators, F11.3/F9.8 and W33, of the sub-arcsecond binary systems with those of other 
Class II objects. The silicate feature properties and the dust processing indicators of objects with 
close (< 1") companions are statistically similar to those of wide binaries and single stars. These 
results suggest that the disk properties are only weakly affected by the presence of companions. 
The multiple systems in our sample have companions with projected separations > 20 AU (see 
Table [3]), and the above results hold for such systems. However, only 12 out of the 68 Class II 
objects in our sample are known to have sub-arcsecond companions. The small number makes the 
statistical comparison less certain. In addition, the multiplicity census of Cha I sample is likely to 
be incomplete, as it is harder to rule out extremely close 1 ") companions. Therefore, we are 
not able to draw any definite conclusions on the effect of multiplicity on the disk structure from the 
comparison based on our sample. Below we discuss the disk properties of Class II objects which 
show clear evidence for disk truncation and examine if close companions are responsible for it. 

We first examine the multiplicity of TDs and PTDs which have holes/gaps within their disks. 
The SED of a close (separation ~ a few AU) binary syste m with a circunibinary disk can appear 
TD-like as has been shown in the case of CoKu Tau/4 (jlreland &: Kraud l2008l ) . However, only 
two of the TDs/PTDs in Cha I, viz., CS Cha and T54 are known to have close stellar companions 
(see Table [3]). T54 has a 0.26" companion (~ 43 AU at 165 pc) and the disk around it could be 
truncated from inside by the gravitational infiuence of the close companion; the mid-IR excess that 
we observe could be from the circumbinary disks surrounding t hese systems. On the o ther hand, 



CS C ha, which has an estimated inner hole of 43 AU i n radius (lEspaillat et al.l l2007al : iKim et al 



20091 ). is a spectroscopic binary (jCuenther et al.l 120071 ) . It is not clear, though, if such a clos e 
companion can carve out the hole of this size in the disk (see lKim et al.ll2009l : iFurlan et al.ll2009l ). 
CHXR 22E, T25, T35, and T56 do not have known companions around them: companions up to 
5,12,12 and 15 times as faint as the primary in the Kg band have been ruled out for sep arations > 
0.1" (~17 AU) for CHXR 22E, T25, T35, and T56 respectively (|Lafreniere et ahlboosl ). Further, 
the compan ions reported in the literature for SZ Cha are not confirmed members of the Cha I 
association (jLuhmanll2007l ) . Comparative studies of various mechanisms that could possibly create 
holes or gaps in disks have shown that most of the observed characteristics of a large sample of 
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TDs/PTDs from various star fo rming regions are best explained by Jovian mass planets opening 



up gaps and holes iri the di sks (jKim et alj l2009l : iNaiita et al.l 120071 : iRice et al.l 120031 : ICalvet et al 
2OO5I : lEspaillat etaPboOTbl l. 



Next, we discuss objects which show clear evidence for disk truncation, i.e, continuum indices 
bluer than that for an optically thick, flat disk, to see if known companions are responsible for disk 
truncation. Among the Class II disks in our sample, 9 show ne_i3 < -1.33, steeper than that for an 
optically thick, flat disk. 5 of them are TDs, only two of which , CS Cha and T54, are sub-arcsecond 
binaries as discussed above. Among the remaining 4 objects, 2M J11062942-7724586, CHXR 30A, 
C7-1 and ISO 282, only CHXR 30A is known to have a close (~ 0.46") companion which could 
be responsible for the disk truncation. Two class II objects, T51 and Hn 5, show ni3_3i < - 
1.33, indicating that they have outwardly truncated d isks: no sub-arcseco nd companions have been 
reported for these two objects (see Section l4?T]l (see iFurlan et al.ll2009l ). Thus, there is evidence 
that disk truncation is caused by a close companion in a few objects; however, no sub-arcsecond 
companions are known for most objects with truncated disks. While undetected companions at 
small separations (^ 1") could still be responsible for disk truncation in these objects, available 
multiplicity data for Class II objects in Cha I do not clearly demonstrate that close (projected 
separation > 20 AU) companions affect disk structure. Similar results have also been found for 
Class II disks in Taurus, where no statistically significant differences were foun d in the structure 
and dust properties of disks in single and multiple systems (jPascucci et al.ll2008l ). 



6. Accretion and disk structure 

In Section 3, we classified the Class II objects in Cha I into WTTS and CTTS using the Hq 
emission strength as an accretion indicator (see Table [2]). WTTS are young T Tauri stars that have 
little if any accretion and probably have dissipated their inner disk, unlike CTTS that are actively 
accreting from the disk. Here we compare the mid-IR emission of the WTTS and CTTS in our 
sample as revealed by their IRS spectra. There are 48 CTTS and 15 WTTs in our Cha I sample: 
Two of the WTTS are Class III objects and two are TDs. 

All of the WTTS in our sample except the Class III objects show mid-IR (> 6 /im) excess 
emission as can be seen from their SEDs (Figure [5]). All of them except the 2 TDs, CHXR 22E 
and T25, and CHSM 10862, also show both 10 and 20 //m silicate feature in emission. The mid-IR 
excess and the silicate features seen in the IRS spectra of the WTTs indicate that they still have 
dusty disks around them. However, most WTTS show low or no near-IR (< 6 /im) excess emission 
(Figure [5]). In Figure \32\ we compare the distribution of spectral indices n2_6, ?t.6-13 and ni3_3i 
for WTTS and CTTS. The distribution of all three continuum indices for the WTTS is shifted to 
more negative values compared to those for the CTTS, indicating that the WTTS have steeper 
continuum slope of the SED in the infrared wavelengths. As can be seen from Figure \32\ this effect 
is more conspicuous at shorter wavelengths, in the n2-Q index than in the other two indices. A two- 
sided K-S test shows that the distribution of the n2-6 index of CTTS and WTTS are statistically 
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different at 99.99% significance level. 



Most TDs and PTDs in Cha I, except CHXR 22E and T25, could be classified as CTTS based 
on the strength of their Ha emission (see Table [2]). T54, for which Ha measurement is not available 
i n the literature, has been reported to be accreting at a relatively high rate based on its U band flux 
(jKim et al.ll2009l ). Thus most TDs/PTDs is Cha I are active accretors. All the objects identified as 
having enhanced 10 /im emission (Wio > ^ lo, upper octUe] see Figure [23]) . except T33A, are CTTS 
based on their Ha equivalent widths indicating that they are actively accreting from their disks. 
T33A is a close binary and the IRS spectra has contribution from both the components. The 
interpretation of Wio and the continuum index in this system is inconclusive. Thus most objects 
in Cha I which show evidence for holes or gaps in their disks appear to be actively accreting. 

Apart from the TDs and PTDs discussed earlier, there are a few other objects in our sample 
which show significant fiux deficits at wavelengths < 8 ^um, but show silicate emission features and 
continuum excesses longward of 8 fim with ni3_3i index comparable to that of Class H objects. 
They are CHXR 20, CHXR 30A, C7-1, and CHXR 47. The slope of their continuum between 2 and 
13 fim is steeper (more negative) than that for most (~ 80%) Class H disks and they show little 
excess emission below 8 /im. The dust grains in their inner disks appear to be greatly depleted 
creating an opacity deficit. They differ from the TDs and PTDs in that their excess flux, even at 
wavelengths longer than 13 /im, is considerably lower than that of the Class II median. Also, they 
do not have a rising continuum longward of 13 fim indicating that the outer disk probably is not as 
optically thick as that of Class II disks and the 'wall' emission is not significant. This class of disks 
has been identifi ed by various names in the literature: 'ev o lved disks', 'anemi c disks', 'homologousl y 



depleted disks' (jLada et al 



2006 



Hernandez et al. 



2007 



Currie et al.ll2009l : iLuhman et al.ll2010l ). 



Among these objects, CHXR 20, CHXR 30A and CHXR 47 are WTTS (see Tabled]); they have 
stopped accreting or are accreting at very low rates. C7-1, however, has an Ha equivalent width 
of 60 A (lLuhmanll2004al ) and is accreting. It probably has gas in the inner most parts of the disk. 
It is not clear if evolved disks fit into an evolutionary sequence in which the optically thick disk 
develops inner gaps and holes and eventually dissipates inside out. It is possible that they represent 
a different evolutionary path for the primordial optically thick disks where the disk dissipation does 
not proceed inside out. In this sense, evolved disks could also be considered to be in a transitional 



phas e between the prini o rdial optically thick d isks and the optically thin debris disks (jLada et al 



2006 : ICurrie et al.l l2009l : iMuzerolle et al.ll20ld ). A more detailed discussion of these objects and 
the question of disk dissipation in Cha I will be addressed in a forthcoming paper (Furlan et al. in 
preparation). 



7. Summary and conclusions 

We have presented the mid-IR spectra of 82 young sources (most of them later than K5, 
including very low mass stars and brown dwarfs) in the Chamaeleon I star forming region, obtained 
with the IRS onboard the Spitzer Space Telescope. The IRS spectra and the photometry compiled 
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from the literature have been used to construct the SEDs of these objects. Based on the observed 
SEDs, we classified the young stars in Cha I into various evolutionary classes. Our sample consists 
of 6 Class I objects, 5 Flat spectrum sources, 68 Class II objects and 3 Class III stars. Eight of the 
Class II sources appear to be transitional disk candidates. 

We analyzed the IRS spectra of the Class II objects in order to study the structure and 
evolution of their disks. To do this, we compared the shape of the observed mid-IR continuum to 
that predicted by the irradiated, accretion disk models which incorporates effects of disk evolution, 
like grain growth and sedimentation. We analyzed the silicate emission features in the IRS spectra 
to study the structure and dust properties of protoplanetary disks in Cha I. We computed median 
SEDs representative of T Tauri stars (K5-M2) and very low mass objects (M3-M8) in our sample 
and compared them to the model SEDs and to each other to investigate the differences in disk 
properties as a function of stellar mass. We also examined the effects of multiplicity and accretion 
on the observed disk properties. Finally, we compared the properties of protoplanetary disks in 
Cha I (age ~ 2 Myr) with those in slightly younger (ages < 1 Myr) Taurus and Ophiuchus star 
forming regions. The main results of this study are summarized below. 



• The observed continuum indices ng-is and ni3_3i for most Class II objects are best explained 
by the disk models where the surface layers of the disk are depleted of small grains by a factor 
of 100 to 1000 compared to that of the standard ISM dust-to-gas mass ratio. Substantial dust 
settling has occurred in most protoplanetary disks surrounding Class II objects in Cha I. The 
distribution of ng-is and ni3_3i continuum indices of Class II objects in the ~ 2 Myr old Cha 
I, ~ 1 Myr Taurus and < 1 Myr Ophiuchus are statistically indistinguishable indicating that 
the degree of dust settling in the disks and the frequency of settled disks are quite similar in 
these three regions. Most protoplanetary disks undergo significant dust settling by < 1 Myr. 
No evidence for temporal evolution is seen in the observed degree of settling in < 1 — 2 Myr 
old disks. Grain growth and sedimentation in protoplanetary disks appear to be balanced by 
counteracting processes such as fragmentation of large grains in the midplane and turbulent 
mixing. 

• The absence of objects showing high n2-6 index in Cha I indicates a paucity of high accretors 
among the Class II objects in Cha I as compared to Taurus and Ophiuchus. The apparent 
lack of high accretors in Cha I is consistent with the accretion rates in Class II disks declining 
with age. 



We confirmed the TDs and PTDs in Cha I reported by iKim et al.l (|2009l ) and identified and 
additional TD candidate of spectral type M5, 2M J11241186- 7630425. TDs and PTDs show 
presence of holes or gaps in their inner disks, but have optically thick outer disks similar to 
that in Class II disks. Most TDs and PTDs are active accretors. Some of them have close 
stellar companions which are likely responsible for their altered disk structure; but in most 
TDs and PTDs, the holes and gaps are probably caused by very low mass, hitherto undetected 
substellar or planetary mass objects. 



- 26 - 



We also identified a few objects which show significantly high 10 fj,m silicate equivalent widths 
(Wio > ^ 10, upper octiie) than what the models for 'full' disks would predict, but have ni3_3i 
index within the al lowable range f or the 'full' disks. These objects are possibly opening up 



gaps in their disk (jEspaillatl l2009l : iFurlan et al.ll2009l ). However, detailed modeling of their 



disk structure is required to confirm this. Most of these objects are actively accreting. 

We further confirmed 4 objects which could be classified as 'evolved' disks based on their 
observed SEDs. Their continuum emission is considerably lower than the 'full ' disks in the 
entire infrared wavelength range (2 - 36 /xm) indicating that their outer disks are relatively 
optically thin. Most of these objects have stopped accreting or are accreting at a low rate. 

• Almost all Class II disks show mid-IR continuum excess and silicate emission features centered 
at 10 and 20 fim. The strength of both the 10 and 20 ^m silicate features, on average, decreases 
as the degree of dust settling increases. Flared disks with less sedimentation show stronger 
silicate features whereas flatter disk with high degree of dust settling show weaker silicate 
features. We further find that the relative strength of the 20 fim feature over the 10 /im 
feature drops with increasing degree of dust sedimentation in the disks. 

• The 10 and 20 /xm silicate features in Class II disks display a wide variety of feature strengths 
and shapes. The flux ratio Fn.s/Fg.g and the equivalent width of the 33.6 fim forsterite feature 
indicate that the dust grains in the disk in Cha I have undergone significant processing (grain 
growth and crystallization). Additionally, we find a tight correlation between the equivalent 
width of 10 fim silicate feature and the flux ratio Fn.s/Fg.g, which is a measure of the degree 
of dust processing in the inner (< 1-2 AU) disks. Less settled disks, with stronger silicate 
emission show low degree of dust processing. The more settled disks in our sample, in general, 
show higher degree of dust processing. 

However, disks with radial holes/gaps in general, show larger silicate emission equivalent 
widths and relatively pristine dust composition compared to that found for radially continuous 
'fuh' disks. 

We find a tight correlation between the flux ratio Fn.a/Fg.g and the equivalent width of the 
forsterite feature at 33.6 /im which suggests that the inner (< 1-2 AU) and outer (< 10 AU) 
disk crystallinity track each other well for Class II disks in Cha I. The correlation between 
crystalline content of the inner and outer regions of the disks has also been found for Class II 
objects in Taurus and Ophiuchus region. 

• Median SEDs computed for Class II objects in the spectral type range of K5-M2 and M3-M8 
show crucial difference in the shape and strength of the continuum emission at wavelengths 
< 8 //m: mid-M type stars and brown dwarfs show excess emission only at wavelengths > 
4-5 iim while in T Tauri stars (K5-M2) the excess emission begins at ~ 2 /im; the observed 
continuum at wavelengths < 8 i-im in mid-M type stars and brown dwarfs is considerably 
steeper than that observed for T Tauri stars. However, both the medians look similar at 
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longer wavelengths: they show 10 and 20 /xm sihcate emission features and their continuum 
slopes at wavelengths > 8 /xm are not very different. 

• Most WTTS in our sample show mid-IR (> 8 /um) excess emission and the silicate features 
at 10 and 20 /xm, indicating the presence of disks around them. However, most of them 
do not show measurable excess at near-IR (< 6 fim) wavelengths. Their mid-IR continuum 
indices are steeper (bluer) than those for CTTS. The observed mid-IR excess of WTTS is 
considerably weaker than that of CTTS at all wavelengths. 
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Table 1. Log of IRS observations of the Cha I sample 



Object Name 2MASS ID RA Dec Date of Campaign AOR Module 

(2000) (2000) observation No. ID 

(1) (2) (3) (4) (5) (6) (7) (8) 



SX Cha 


J10555973-7724399 


10 


55 59.73 


-77 


24 39.90 


2005 


May 26 


21 


12697345 


SLSHLH 


T5 


J10574219- 7659356 


10 


57 42.20 


-76 


56 35.70 


2005 


Jul 13 


22 


12696320 


SLLL 




J10580597- 7711501 


10 


58 5.98 


-77 


11 50.10 


2007 


Jun 26 


41 


21914368 


SLLL 


SZ Cha 


J10581677- 7717170 


10 


58 16.77 


-77 


17 17.10 


2005 


Jul 12 


22 


12696832 


SLLL 


TW Cha 


J10590108- 7722407 


10 


59 1.09 


-77 


22 40.70 


2005 


Jul 10 


22 


12696576 


SLLL 


CR Cha 


J10590699- 7701404 


10 


59 6.99 


-77 


1 40.40 


2005 


May 26 


21 


12697345 


SLSHLH 




J11011926-7732383 


11 


1 


19.22 


-77 


32 38.60 


2007 


Jun 18 


41 


21920512 


SLLL 


CS Cha 


J11022491-7733357 


11 


2 


24.91 


-77 


33 35.70 


2005 


Jul 11 


22 


12695808 


SLLL 


CT Cha 


J11040909- 7627193 


11 


4 


9.09 


-76 


27 19.40 


2005 


May 26 


21 


12697345 


SLSHLH 


T14a 


J11042275- 7718080 


11 


4 


22.76 


-77 


18 8.10 


2005 


Jul 11 


22 


12695808 


SLLL 


ISO 52 


J11044258- 7741571 


11 


4 


42.58 


-77 


41 57.10 


2005 


Jul 12 


22 


12691200 


SLLL 


T21 


J11061540- 7721567 


11 


6 


15.41 


-77 


21 56.80 


2005 


Jul 13 


22 


12696320 


SLLL 




J11062942- 7724586 


11 


6 


29.43 


-77 


24 58.60 


2006 


Aug 2 


33 


18360832 


SLLL 


CHSM 7869 


J11063276-7625210 


11 


6 


32.77 


-76 


25 21.10 


2007 


Jun 26 


41 


21913600 


SLLL 


ISO 79 


J11063945-7736052 


11 


6 


39.45 


-77 


36 5.20 


2005 


Jul 11 


22 


12690176 


SLLL 


Hn 5 


J11064180-7635489 


11 


6 


41.81 


-76 


35 49.00 


2005 


Jul 13 


22 


12696320 


SLLL 


UX Cha 


J11064346- 7726343 


11 


6 


43.47 


-77 


26 34.40 


2005 


Jul 13 


22 


12696064 


SLLL 


CHXR 20 


J11064510- 7727023 


11 


6 


45.10 


-77 


27 2.30 


2005 


Jul 12 


22 


12695296 


SLLL 


CcdllO IRS4 


J11064658-7722325 


11 


6 


46.58 


-77 


22 32.60 


2005 


Apr 17 


20 


12692224 


SLLL 


ISO 86 


J11065803-7722488 


11 


6 


58.03 


-77 


22 48.90 


2005 


Apr 17 


20 


12692224 


SLLL 


UY Cha 


J11065906-7718535 


11 


6 


59.07 


-77 


18 53.60 


2005 


Apr 17 


20 


12692224 


SLLL 




J11065939-7530559 


11 


6 


59.40 


-75 


30 56.00 


2007 


Jun 15 


41 


21915136 


SLLL 




J11070369- 7724307 


11 


7 


3.69 


-77 


24 30.70 


2006 


Aug 2 


33 


18360576 


SLLL 


CedllO IRS6 


J11070919- 7723049 


11 


7 


9.20 


-77 


23 5.00 


2005 


Jul 12 


22 


12696832 


SLLL 


ISO 91 


J11070925- 7718471 


11 


7 


9.25 


-77 


18 47.10 


2005 


Apr 24 


20 


12691712 


SLLL 


UZ Cha 


J11071206-7632232 


11 


7 


12.07 


-76 


32 23.20 


2005 


Jul 13 


22 


12696320 


SLLL 


CHXR 22E 


J11071330-7743498 


11 


7 


13.33 


-77 


43 49.80 


2006 


Jul 3 


32 


18361344 


SLLL 


Cha Hal 


J11071668-7735532 


11 


7 


16.69 


-77 


35 53.30 


2005 


Jul 13 


22 


12689920 


SLLL 


T25 


J11071915-7603048 


11 


7 


19.15 


-76 


3 4.80 


2005 


Apr 24 


20 


12695552 


SLLL 


DI Cha 


J11072074-7738073 


11 


7 


20.74 


-77 


38 7.40 


2005 


May 26 


21 


12697345 


SLSHLH 


B35 


J11072142- 7722117 


11 


7 


21.43 


-77 


22 11.80 


2005 


Jul 13 


22 


12696320 


SLLL 


VV Cha 


J11072825-7652118 


11 


7 


28.26 


-76 


52 11.90 


2005 


Jul 13 


22 


12696064 


SLLL 



Table 1 — Continued 



Object Name 2MASS ID RA Dec Date of Campaign AOR Module 

(2000) (2000) observation No. ID 

(1) (2) (3) (4) (5) (6) (7) (8) 





T1 1 07424 '1-7733593 


11 


7 42.45 


-77 


33 


59.40 


2005 


Jul 12 


22 


1 2690688 


SLLL 


T28 


T1 1 074366-773941 1 


11 


7 43.66 


-77 


39 


41.10 


2005 


Jul 12 


22 


12696832 


SLLL 




Tl 1 074656- 7fil 51 74 

O i- -LyJ 1 '-t\J'J\J 1 \J A.IJ 1. 1 '± 


11 


7 46.56 


-76 


15 


17.50 


2007 




41 


21 91 3856 


SLLL 


CHXR SOB 


Tl 1 075730-771 7262 


11 


7 57.31 


-77 


17 


26.20 


2005 


Jul 10 


22 


1 2696576 


SLLL 


T29 


J11075792-7738449 


11 


7 57.93 


-77 


38 


44.90 


2005 


May 26 


21 


12697345 


SLSHLH 


CHXR 30A 


1 11080002-771 7304 


11 


8 0.02 


-77 


17 


30.50 


2005 


Mav 30 


21 


12692481 


SLLL 


VW Cha 


J11080148-7742288 


11 


8 1.49 


-77 


42 


28.80 


2005 


Jul 12 


22 


12696832 


SLLL 


ISO 12fi 


71 1 080297-7738425 


11 


8 2.98 


-77 


38 


42.60 


2005 


1V± <Xy 


21 


1 2697345 


SLSHLH 


CTI Cha 


Tl 1 080329-77391 74 


11 


8 3.30 


-77 


39 


17.40 


2006 


Mar 7 


29 


1 2697088 


ST,SHT,H 


T33A 


Til 081 509-7733531 


11 


8 15.10 


-77 


33 


53.20 


2006 


Mar 7 


29 


1 2697088 


SLSHLH 


ISO 138 


Til 081 850-7730408 


11 


8 18.50 


-77 


30 


40.80 


2005 


Jul 12 


22 


1 2690432 


SLLL 


ISO 143 


Til 082238-7730277 


11 


8 22.38 


-77 


30 


27.70 


2005 


Jul 12 


22 


12690432 


SLLL 


Cha IRN 


J11083896-7743513 


11 


8 38.20 


-77 


43 


51.80 


2005 


May 26 


21 


12697345 


SLSHLH 


T35 


Tl 1083905-7716042 


11 


8 39 05 


-77 


16 


4.20 


2005 


Jul 13 


22 


12696320 


SLLL 


VY Cha 


J11085464-7702129 


11 


8 54.64 


-77 


2 13.00 


2005 


May 26 


21 


12691969 


SLLL 


Cl-6 


J11092266-7634320 


11 


9 22.67 


-76 


34 


32.00 


2005 


May 26 


21 


12697345 


SLSHLH 


VZ Cha 


J11092379-7623207 


11 


9 23.79 


-76 


23 


20.80 


2005 


Jul 12 


22 


12696832 


SLLL 


ISO 192 


J11092855-7633281 


11 


9 28.55 


-76 


33 


28.10 


2006 


Mar 7 


29 


12697088 


SLSHLH 


Cl-25 


J11094192-7634584 


11 


9 41.93 


-76 


34 


58.40 


2005 


Jul 12 


22 


12696832 


SLLL 


C7-1 


J11094260-7725578 


11 


9 42.60 


-77 


25 


57.90 


2005 


Jul 12 


22 


12686336 


SLLL 


Hn lOE 


J11094621-7634463 


11 


9 46.21 


-76 


34 


46.40 


2005 


Jul 13 


22 


12696320 


SLLL 


B43 


J11094742-7726290 


11 


9 47.42 


-77 


26 


29.10 


2005 


Apr 17 


20 


12691456 


SLLL 


HD 97300 


J11095003-7636476 


11 


9 50.03 


-76 


36 


47.70 


2006 


Mar 7 


29 


12697088 


SLSHLH 


ISO 220 


J11095336-7728365 


11 


9 53.37 


-77 


28 


36.60 


2006 


Mar 7 


29 


12686080 


SLLL 


T42 


J11095340-7634255 


11 


9 53.41 


-76 


34 


25.50 


2006 


Mar 7 


29 


12697088 


SLSHLH 


T43 


J11095407-7629253 


11 


9 54.08 


-76 


29 


25.30 


2005 


Apr 17 


20 


12692224 


SLLL 


ISO 225 


J11095437-7631113 


11 


9 54.38 


-76 


31 


11.40 


2005 


Jul 12 


22 


12690432 


SLLL 


Cl-2 


J11095505-7632409 


11 


9 55.06 


-76 


32 


41.00 


2005 


Jul 12 


22 


12696832 


SLLL 


WX Cha 


J11095873-7737088 


11 


9 58.74 


-77 


37 


8.90 


2005 


Jul 10 


22 


12696576 


SLLL 


WW Cha 


JlllOOOlO-7634578 


11 


10 0.11 


-76 


34 


57.90 


2006 


Mar 7 


29 


12697088 


SLSHLH 


Hn 11 


J11100369-7633291 


11 


10 3.69 


-76 


33 


29.20 


2005 


Jul 13 


22 


12696064 


SLLL 


T45a 


J11100469-7635452 


11 


10 4.69 


-76 


35 


45.30 


2005 


May 26 


21 


12691969 


SLLL 



Table 1 — Continued 



Object Name 2MASS ID RA Dec Date of Campaign AOR Module 

(2000) (2000) observation No. ID 



(1) 


(2) 




(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


WY Cha 


J11100704- 7629376 


11 


10 


7.04 


-76 


29 37.70 


2005 Jul 13 


22 


12696064 


SLLL 


ISO 235 


Jl 1100785-7727480 


11 


10 


7.85 


-77 


27 48.10 


2005 Jul 12 


22 


12690432 


SLLL 


ibO 26" 


J 11101141-7635292 


11 


10 


11.42 


-76 


35 29.30 


2005 Jul 10 


22 


12696576 


SLLL 


/"ITT'VT^ A '-7 

OHXR 47 


J 1 1 103801-7732399 


11 


10 


38.02 


-77 


no ctr\ 

32 39.90 


2005 Jul 13 


22 


12696064 


SLLL 


ISO 252 


Jl 1104141-7720480 


11 


iU 




-77 


20 48.00 


2005 Jul 9 


22 


12685825 


SLLL 


T47 


J11104959- 7717517 


11 


10 


49.60 


-77 


17 51.70 


2005 Jul 10 


22 


12696576 


SLLL 


WZ Cha 


J11105333-7634319 


11 


10 


53.34 


-76 


34 32.00 


2005 Apr 17 


20 


12692224 


SLLL 


ISO 256 


J11105359- 7725004 


11 


10 


53.59 


-77 


25 0.50 


2005 May 26 


21 


12691969 


SLLL 


Hn 13 


.111105597-7645325 


11 


10 


55.97 


-76 


45 32.60 


2005 Jul 12 


22 


12690688 


SLLL 


XX Cha 


J11113965-7620152 


11 


11 


39.66 


-76 


20 15.20 


2005 Jul 10 


22 


12696576 


SLLL 


ISO 282 


J11120,351-7726009 


11 


12 


3.51 


-77 


26 1.00 


2005 Jul 13 


22 


12686848 


SLLL 


T50 


J11120984-7634366 


11 


12 


9.85 


-76 


34 36.50 


2005 Jul 13 


22 


12696064 


SLLL 


T51 


J11122441-7637064 


11 


12 


24.41 


-76 


37 6.40 


2005 Jul 10 


22 


12696576 


SLLL 


CV Cha 


J11122772-7644223 


11 


12 


27.72 


-76 


44 22.30 


2006 Mar 7 


29 


12697088 


SLSHLH 


T54 


J11124268- 7722230 


11 


12 


42.69 


-77 


22 23.00 


2005 Apr 24 


20 


12695552 


SLLL 


Hn 21W 


J11142454-773,3062 


11 


14 


24.54 


-77 


33 6.20 


2007 Jun 16 


41 


21914112 


SLLL 


T56 


J11173700-7704381 


11 


17 


37.01 


-77 


4 38.10 


2005 Jul 13 


21 


12696064 


SLLL 




J11183572- 7935548 


11 


18 


35.73 


-79 


35 54.90 


2007 Jun 10 


41 


21915392 


SLLL 




J11241186-7630425 


11 


24 


11.87 


-76 


30 42.60 


2007 Aug 2 


42 


21914624 


SLLL 




J 11432669-7804454 


11 


43 


26.70 


-78 


4 45.50 


2007 Jun 17 


41 


21914880 


SLLL 



Note. — Columns (1) & (2) list the object names and their 2MASS identifiers. Observed coordinates are listed in Columns (3) 
and (4) and the date of observations in column (5). The Spitzer IRS observing campaigns during which the objects were observed 
are listed in Column (6) and the corresponding Astronomical Observation Requests (AORs) in Column (7). Column (8) lists the 
IRS modules used in the observations. 
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Table 2. Basic data for the Cha I sample 



Object Name 


Sp. Type 


Aj 


"2-25 


"2-26 


SED Class 


"5-12 


W(H„) 


c/w 






(mag) 


(observed) 


(dereddened) 






(A) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SX Cha 


MO 


0.7 


-0.61 


-0.59 


II 


-0.71 


31.7 


C 


T5 


M3.25 


0.3 


-1.02 


-1.02 


II 


-0.87 


22.0 


c 


2M J10580597-7711501 


M5.25 


0.4 


-0.94 


-0.93 


II 


-1.18 






SZ Cha 


KO 


0.4 


-0.59 


-0.58 


II/T 


-1.20 


7.0 


c 


TW Cha 


K8 


0.3 


-1.00 


-1.00 


II 


-1.35 


34.5 


c 


CR Cha 


K2 


0.0 


-0.82 


-0.82 


II 


-0.62 


38.5 


c 


2M J11011926- 7732383 


M7.25 


0.4 


-2.38 


-2.37 


III 


-2.70 


25.0 


w 


CS Cha 


K6 


0.0 


-0.64 


-0.63 


II/T 


-1.52 


20.0 


c 


CT Cha 


K5 


0.4 


-0.64 


-0.63 


II 


-0.50 


34.5 


c 


T14a 


K7 


0.4 


0.50 


0.51 


^ 


0.52 


159.0 


c 


ISO 52 


M4 


0.3 


-0.81 


-0.81 


II 


-0.47 


6.4 


w 


T21 


G5 


0.9 


-2.53 


-2.50 


III 


-2.65 






2M J11062942- 7724586 


M6 


5.6 


-1.36 


-1.06 


II 


-1.29 






CHSM 7869 


M6 


0.4 


-1.09 


-1.08 


II 


-1.13 


120.0 


c 


ISO 79 


M5.25 


0.7 


-0.70 


-0.69 


II 


-0.76 


8.0 


w 


Hn 5 


M4.5 


0.3 


-0.83 


-0.82 


II 


-0.72 


48.0 


c 


UX Cha 


M3 


0.9 


-2.62 


-2.54 


III 


-2.56 


3.0 


w 


CHXR 20 


K6 


1.0 


-1.30 


-1.27 


II 


-0.81 


1.8 


w 


CedllO IRS4 






1.61 






0.31 






ISO 86 






0.91 






-0.95 






UY Cha 


M4.25 


0.0 


-1.00 


-1.00 


II 


-0.52 


8.5 


w 


2M J11065939- 7530559 


M5.25 


0.1 


-1.00 


-0.99 


II 


-1.23 


51.0 


c 


2M J11070369- 7724307 


M7.5 


4.6 


-0.94 


-0.70 


II 


-0.75 






CedllO IRS6 






0.70 




^ 


-0.81 






ISO 91 


M3 


4.5 


-1.15 


-0.91 


II 


-0.53 






UZ Cha 


MO. 5 


0.6 


-1.04 


-1.03 


II 


-1.02 


9.7 


w 


CHXR 22E 


M3.5 


1.3 


-1.79 


-1.74 


II/T 


-2.62 


4.5 


w 


Cha Hal 


M7.75 


0.0 


-0.52 


-0.52 


II 


-0.18 


112.0 


c 


T25 


M2.5 


0.4 


-0.81 


-0.80 


II/T 


-1.97 


5.0 


w 


DI Cha 


G2 


0.7 


-0.99 


-0.98 


II 


-1.12 


18.3 


c 


B35 






-0.02 




FS 


-0.02 






VV Cha 


M3 


0.1 


-0.94 


-0.94 


II 


-0.40 


72.0 


c 


Cha Ka2 


M5.25 


1.0 


-1.28 


-1.24 


II 


-1.07 


60.0 


c 


T28 


MO 


1.3 


-1.11 


-1.07 


II 


-1.01 


78.0 


c 


CHSM 10862 


M5.75 


0.4 


-0.84 


-0.83 


II 


-0.90 


18.0 


w 


CHXR SOB 


M1.25 


3.1 


-0.96 


-0.79 


II 


-1.42 


80.0 


c 


T29 


K6 


1.2 


-0.50 


-0.47 


II 


-0.36 


61.5 


c 


CHXR 30A 


K8 


2.9 


-1.75 


-1.59 


II 


-1.58 


1.0 


w 


VW Cha 


K8 


0.7 


-0.89 


-0.87 


II 


-1.02 


64.0 


c 


ISO 126 


M1.25 


1.3 


-0.13 


-0.09 


FS 


-0.49 


120.0 


c 


CU Cha 


B9.5 


0.7 


-0.14 


-0.13 


II 


-0.03 






T33A 


G7 


0.8 






II 


-0.47 


0.9 


w 


ISO 138 


M6.5 


0.0 


-1.00 


-1.00 


II 


-0.75 


19.0 


w 
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Table 2 — Continued 



Object Name 


Sp. Type 


Aj 


"2-25 


"2-25 


SED Class 


"5-12 


W(Hc) 


c/w 






(mag) 


(observed) 


(dereddened) 






(A) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ISO 143 


M5 


1.1 


-1.16 


-1.12 


II 


-0.85 


118.0 


C 


Cha IRN 






1.24 




I 


1.46 






T35 


K8 


1.3 


-1.23 


-1.19 


II/T 


-2.48 


100.0 


c 


VY Cha 


M0.5 


0.9 


-0.91 


-0.88 


II 


-1.18 


55.0 


c 


Cl-6 


M1.25 


3.2 


-0.64 


-0.46 


II 


-0.50 


110.0 


c 


VZ Cha 


K6 


0.5 


-1.07 


-1.06 


II 


-0.94 


41.3 


c 


ISO 192 






1.24 




I 


1.49 






Cl-25 






-0.26 




FS 


-0.11 






C7-1 


M5 


1.9 


-1.59 


-1.53 


II 


-1.78 


60.0 


c 


Hn lOE 


M3.25 


1.0 


-0.44 


-0.41 


II 


-0.14 


68.0 


c 


B43 


M3.25 


2.2 


-0.97 


-0.90 


II 


-1.08 


200.0 


c 


HD 97300 


B9 


0.7 


-1.32 


-1.31 


II 


0.17 






ISO 220 


M5.75 


1.7 


-1.23 


-1.18 


II 


-0.94 


150.0 


c 


T42 


K5 


1.4 


-0.39 


-0.34 


II 


-0.38 


205.0 


c 


T43 


M2 


1.4 


-0.94 


-0.89 


II 


-0.92 


43.0 


c 


ISO 225 


Ml. 75 


1.2 


0.10 


0.14 


FS 


0.00 


12.7 


c 


Cl-2 






0.10 




FS 


-0.43 


45.0 




WX Cha 


Ml. 25 


0.5 


-1.06 


-1.05 


II 


-1.16 


63.0 


c 


WW Cha 


K5 


1.3 


-0.46 


-0.42 


II 


-0.78 


61.0 


c 


Hn 11 


K8 


2.1 


-0.84 


-0.77 


II 


-0.58 


19.0 


c 


T45a 


Ml 


0.5 


-1.00 


-1.00 


II 


-0.93 


7.5 


W 


WY Cha 


MO 


1.1 


-1.18 


-1.15 


II 


-0.83 


51.0 


c 


ISO 235 


M5.5 


2.4 


-1.31 


-1.24 


II 


-0.71 


18.0 


W 


ISO 237 


K5.5 


1.9 


-0.78 


-0.72 


II 


-0.88 


3.8 


c 


CHXR 47 


K3 


1.4 


-1.24 


-1.20 


II 


-1.05 


1.0 


W 


ISO 252 


M6 


0.9 


-1.01 


-0.98 


II 


-0.71 


320.0 


c 


T47 


M2 


1.1 


-0.63 


-0.60 


II 


-0.19 


89.5 


c 


vv LI v^iia 


iViO. / u 


n 1 


n 7Q 
-u. / y 


n 7Q 
-u. / y 


jj 


-u.oo 


97Q n 
z / y.u 


n 


ISO 256 


M4.5 


2.5 


-0.75 


-0.67 


II 


-0.57 


55.0 


c 


Hn 13 


M5.75 


0.2 


-1.02 


-1.02 


II 


-1.08 


22.5 


c 


XX Cha 


M2 


0.3 


-0.90 


-0.89 


II 


-0.90 


90.0 


c 


ISO 282 


M4.75 


1.0 


-1.13 


-1.10 


II 


-1.51 


50.0 


c 


T50 


M5 


0.1 


-1.08 


-1.08 


II 


-1.01 


19.5 


W 


T51 


K3.5 


0.0 


-1.08 


-1.08 


II 


-1.04 


2.5 


W 


CV Cha 


G9 


0.4 


-0.78 


-0.77 


II 


-0.79 


62.5 


c 


T54 


G8 


0.6 


-1.68 


-1.67 


II/T 


-2.32 






Hn 21W 


M4 


0.7 


-1.34 


-1.33 


II 


-1.26 


68.0 


c 


T56 


M0.5 


0.1 


-0.93 


-0.93 


II/T 


-1.22 


56.0 


c 


2M J11183572- 7935548 


M4.75 


0.0 


-0.66 


-0.66 


II/T 


-0.67 


11.0 


w 


2M J11241186-7630425 


M5 


0.7 


-0.91 


-0.89 


II/T 


-1.13 


30.0 


c 


2M J11432669-7804454 


M5 


0.2 


-1.48 


-1.47 


II 


-1.77 


75.0 


c 



Note. — Spectral types and extinction values (Aj) listed in Columns (2) & (3) are from l|Luhmanll2007l . 120081 1. 'T' in Co lumn (5) 
denotes transitional disk candidates (see Section 3 .2). Column ( 6) lists Ha eq uivalent widths compiled from the literature (^ Luhmanl 
[2004a; Herbig fc Bell. 1988 : .Comeron et al... 19991 : 1 Waited 1 19921 : ICovino et al]ll997. : .Feigelson k. Kriss 1989: Guenther et al..l2007l 'l 



Positive values indicate Ha in emission. 



Table 3. Objects in multiple systems 



Separation Brightness 
Name N ( ) Spectral Type Comparison References 

(1) (2) (3) (4) (5) (6) 



SX Cha (A,B) 


2 


2.2 


M0+M3.5 


A>B 


1,2 


T5 (A,B) 


2 


0.16 


M3.25 


A^B 


2 


2M J11011926-7732383 (A,B) 


2 


1.4 


M7.25+M8.25 


A>B 


3 


Una (A.B) 


2 


SB 


K6 


A>B 


4 


T14a (A,B) 


2 


2.35 


K7 


A>B 


5 


T21 (A,B) 


2 


0.14 


G5 


A>>B 


2 


CHXR 20 + UX Cha 


2 


28.46 


K6+M3 


CHXR20>UX Cha 


6 


CedllO IRS6 (A,B) 


2 


1.95 


ClassI 


A>B 


5,7 


DI Cha (A,Ba,Bb) 


3 


4.68(A,Ba+Bb),0.07(Ba,Bb) 


G2(A),M4.25(Ba+Bb) 


A>>Ba+Bb,Ba~Bb 


1,8,2 


VV Cha (A,B) 


2 


0.8 


M3 


A>B 


1,2 


L^ria rlaz [A^H) 


2 


0.17 


M5.25 


A^B 


9,2 


T28 (A,B) 


2 


28.87 


M0+M5.75 


A>>B 


6 


T29 + ISO 126 (A,B) 


3 


17.6, 0.29(ISO 126 A,B) 


K6+M1.25 


T29>ISO 126A>ISO 126B 


8 


CHXR 30 (Aa,Ab,B) 


3 


0.46(Aa,Ab),9.93(A,B) 


K8(Aa+Ab)+Ml. 25(B) 


Aa>B,Aa>>Ab 


6,2 


VW Cha (A,Ba,Bb,C,D) 


4 


0.66(A,Ba),0.7(A,Bb),0.1(Ba,Bb),2.7(A,C),16.9(A,D) 


K8(A+B),K7(Ba+Bb),M2.5(D) 


A>B,Ba>Bb,A>C,A>>D 


1,8,10,11,12,2 


T33 (A,B) 


2 


2.48 


G7+K6 


A>B 


10,11,5,2 


ISO 143 + ISO 138 


2 


18.16 


M5+M6.5 


ISO 143>ISO 138 


6 


Hn lOE + Cl-25 


2 


19.17 


M3.5 


Hn 10E~Cl-25 


6,13 


HD 97300 (A,B) 


2 


0.82 


B9 


A>>B 


8,2 


T43 (A,B) 


2 


0.8 


M2 


A>B 


2 


WX Cha (A,B) 


2 


0.78 


Ml. 25 


A>>B 


8,2 


WY Cha (A,B) 


2 


0.12 


MO 


A>B 


2 


ISO 237 + T45a 


2 


28.32 


K5.5+M1 


ISO 237>T45a 


6 


CHXR 47 (A.Bj^- 


2 


1.2(2), 0.18(13) 


K3 


A~B(2),A>B(13) 


8,2 


T47 (A,B) 


2 


12.09 


M2+M4.5 


A>B 


6 


Hn 13 (A,B) 


2 


0.13 


M5.75 


A~B 


9,2 


XX Cha (B,A) 


2 


24.38 


M2+K6 


B<A 


6 


T51 (A,B) 


3 


1.85 


K3.5 


A>>B 


1,8,2 


CV Cha + CW Cha 


2 


11.55 


G9+M1 


CV Cha>>CW Cha 


1,8 


T54 (A,B) 


2 


0.26 


G8 


A>>B 


8,2 


Hn 21 (W,B) 


2 


5.52 


M4+M5.75 


W>E 


10,6,2 



Note. — When the companion(s) have common names, we list those names in Column(l). When there are more than one measurement for the separations between 
components, and if they are consistent, we have listed the mean value in CoIumn(3) and cited all the references. The spectral types of the primary and of the other 
components listed in Column (4) are from [Luhma 

3 120071 ). except for VW Cha and 2M J11011926- 7732383, for which they are from the corresponding references listed. 
The relative brightness of the components listed in Column (6) is based on the K-band flux. 

^ The separation and the position angle of the CHXR 47 binary system reported by lGhez et all il997|) and |Lafreniere et all i2008|) arc very different. Wc have separately 
listed the binary separation and flux ratio from both the authors with the references in parenthesis. 



Re fere nces. — (DlReipurth &z ZinneckeHll993l : ('2) iLafreniere et al.1|200a: fS) |Luhman||2004d : (4) |Guenther et al.||2007l: f 5) iHaisch et al.||2004|; f6)|Kraus &: Hillenbrand 
2001 : ('7'> |Persi et al.|200ll : ("S'j iGhez et al.|199l : Q'i lAhmic et al.|200l : qo lBrandner et al.|1996l : flD lBrandeker et al.|20ai| : a2'} |Brandner Zinneckei|199l ; (13) iHartiganl 



1995 
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Table 4. Continuum spectral indices of Class II and Class III objects 



Ol^lPpt TIPITTIP 


71 o 
112' 


-6 




-6 




13 


One- 


13 ' 




31 <^ni3- 


31 


SX Cha 


-0, 


.81 


0. 


.05 


-0, 


.65 


0. 


,05 


-0, 


,49 


0. 


,07 


T5 


-1, 


.73 


0, 


.06 


-0, 


.87 


0. 


,08 


-0, 


,32 


0. 


,07 


2M J10580597-7711501 


-2 


.08 


0, 


.08 


-0 


.93 


0. 


,11 


0, 


,26 


0, 


,09 


SZ Cha 


-2 


.07 


0, 


,08 


-1, 


.04 


0. 


,09 


1, 


,72 


0. 


,03 


TW Cha 


-1 


.56 


0, 


,08 


-1, 


.27 


0. 


,10 


-0, 


,17 


0. 


,10 


CR Cha 


-1 


.82 


0, 


,06 


-0, 


.48 


0. 


,06 


-0, 


,12 


0. 


,04 


2M J11011926-7732383 


-2, 


.11 


0. 


,08 


-2, 


.79 


0. 


,29 


-1, 


,50 


0. 


,81 


CS Cha 


-2, 


.62 


0, 


,07 


-1, 


.46 


0. 


,13 


2, 


,90 


0. 


,11 


CT Cha 


-1 


.15 


0, 


,03 


-0, 


.46 


0, 


,02 


-0, 


,32 


0. 


,07 


ISO 52 


-1 


.98 


0, 


,07 


-0 


,42 


0. 


,08 


0, 


,15 


0. 


,05 


T21 


-2, 


.70 


0, 


,09 


-2, 


,71 


0. 


,14 


-1, 


,07 


0. 


,14 


2M J 11062942-7724586 


-1 


.61 


0, 


,11 


-1, 


.40 


0. 


,27 


-0, 


,16 


0. 


,27 


CHSM 7869 


-1 


.81 


0, 


,07 


-1, 


.01 


0. 


,13 


-0, 


,56 


0. 


,22 


ISO 79 


-1 


.21 


0. 


,08 


-0, 


.71 


0. 


,11 


-0, 


,04 


0. 


,06 


Hn 5 


-0 


.67 


0, 


,05 


-0, 


.77 


0. 


,08 


-1, 


,44 


0. 


,12 


UX Cha 


-2 


.64 


0, 


,10 


-2, 


.66 


0. 


,44 


-2, 


,01 


0. 


,94 


CHXR 20 


-2, 


.15 


0. 


,07 


-0, 


.78 


0. 


,08 


-0, 


,96 


0. 


,06 


UY Cha 


-1 


.70 


0. 


,07 


-0, 


.48 


0. 


,09 


-0, 


,92 


0. 


,07 


2M J11065939-7530559 


-1 


.75 


0, 


,09 


-1, 


.27 


0. 


,16 


0, 


,37 


0. 


,19 


2M J 11070369-7724307 


-1 


.55 


0, 


,12 


-0, 


.68 


0, 


,12 


0, 


,15 


0. 


,07 


ISO 91 


-2 


.38 


0, 


,09 


-0, 


,22 


0. 


,12 


-0, 


,40 


0. 


,11 


UZ Cha 


-1 


.76 


0, 


,08 


-0, 


,95 


0. 


,10 


-0, 


,22 


0. 


,07 


CHXR 22E 


-2, 


.51 


0, 


,09 


-2, 


,70 


0. 


,17 


0, 


,58 


0. 


,14 


Cha Hal 


-1 


.50 


0. 


,08 


0, 


.07 


0. 


,12 


0, 


,01 


0. 


,09 


T25 


-2, 


.75 


0. 


,12 


-1, 


.75 


0. 


,26 


2, 


,80 


0. 


,23 


DI Cha 


-1 


.35 


0, 


,06 


-1, 


.01 


0. 


,06 


-0, 


,67 


0. 


,05 


VV Cha 


-1 


.74 


0, 


,06 


-0, 


.38 


0. 


,07 


-0, 


,78 


0. 


,07 


Cha Ha2 


-1 


.66 


0, 


,08 


-1, 


.04 


0, 


,10 


-0, 


,99 


0. 


,07 


T28 


-1 


.61 


0, 


,06 


-0, 


.99 


0. 


,08 


-0, 


,60 


0. 


,07 


CHSM 10862 


-1 


.41 


0, 


,05 


-0, 


.87 


0. 


,11 


-0, 


,05 


0. 


,11 


CHXR 30B 


-0, 


.81 


0. 


,06 


-1, 


.33 


0. 


,08 


-0, 


,39 


0. 


,09 


T29 


-1, 


.06 


0. 


,03 


-0, 


.07 


0. 


,03 


-0, 


,32 


0. 


,03 



Table 4 — Continued 



Object name n2-6 cr„2_g ne-is (^m-w ^t-is-si CTnis-si 



/^xj"V"T> on A 
UjiAK oUA 


-z.4y 


n no 


-l.ol 


nil 
U.ll 


n ccn 

-U.oy 


n nn 

u.uy 


V VV Una 


1 AQ 


n 

U.Uo 


n no 


U.Uo 


n 1 T 
-U.17 


U.Uo 


KjKJ una 


1 O^J 
-i.ZD 


n ncc 
U.Uo 


n no 
U.Uz 


U.Ud 


1 no 
l.Uo 


n nix 
U.Uo 


1 ooA 






n /I n 
-U.4U 


U.Ub 


-U.oo 


n no 
U.Uo 




-i.by 


0.5o 


-0.69 


n Tr' 

0.76 


n 0*7 

-0.27 


n /I o 

U.4z 




-l.oU 


0.06 


-0.85 


n no 
0.08 


1 n /I 
-1.04 


n 1 n 
U.IU 




-Z.Ul 


c\ no 
0.08 


-Z.56 


n on 
O.zO 


1 /IT 
1.47 


n 1 n 

u.iy 


"\ T'\/' r^\^ « 
V Y Una 




0.06 


-1.13 


n nT 

0.07 


n 1 n 

-0.19 


n nrr 

0.05 


Ul-0 


-U.oy 


0.04 


-0.45 


n n /I 

0.04 


n 1 

-0.61 


n nc 

0.05 


VZ Una 


1 A'X 

-i.4i 


0.06 


nor: 

-U.oo 


n n 

0.05 


n ni 

-0.91 


n n*? 

0.07 


U /-i 


o no 
-z.Uz 


n no 
O.Oo 


1 1 

-i.oi 


n 1 n 
0.10 


n 

-0.56 


nil 
0.11 


Hn iUhj 


1 A'X 

-1.4i 


0.07 


n ncc 
-U.Uo 


n no 
0.08 


n OT 
0.37 


n ncc 
0.05 


"D /I O 

r>4o 


-l.o4 


0.06 


-1.11 


n nn 
0.09 


0.43 


n 1 n 
0.10 


riu y7oUU 


-Z.5U 


c\ c\c\ 

0.09 


-0.15 


n OO 


n o c 

-0.85 


n OO 

0.23 


IdU zzU 


-1.67 


0.07 


n no 

-0.98 


n nn 

0.09 


n OO 

-0.82 


n 1 n 

0.10 


i4z 


-U.YO 


0.04 


not; 

-0.35 


n no 

0.03 


n nn 

0.00 


n no 

0.02 




-1.36 


0.05 


n on 

-0.89 


n n/^ 

0.06 


n o'v 

-0.37 


n n r" 

0.05 


W A Cna 


-1.15 


0.06 


-1.14 


n A'? 

0.07 


1 no 

-1.02 


n nT 

0.07 


WW Cha 


-0.90 


0.07 


-0.69 


0.05 


0.42 


0.02 


Hn 11 


-1.05 


0.04 


-0.65 


0.05 


-0.73 


0.05 


T45a 


-2.06 


0.05 


-0.89 


0.07 


0.14 


0.06 


WY Cha 


-1.53 


0.04 


-0.88 


0.06 


-1.16 


0.08 


ISO 235 


-1.69 


0.06 


-0.75 


0.09 


-1.26 


0.09 


ISO 237 


-1.65 


0.06 


-0.81 


0.08 


0.52 


0.06 


CHXR 47 


-1.87 


0.07 


-1.08 


0.08 


-0.66 


0.06 


ISO 252 


-1.40 


0.05 


-0.72 


0.07 


-0.87 


0.09 


T47 


-1.46 


0.04 


-0.04 


0.04 


-0.27 


0.05 


WZ Cha 


-1.29 


0.05 


-0.39 


0.08 


-0.75 


0.08 


ISO 256 


-0.92 


0.07 


-0.65 


0.08 


-0.62 


0.07 


Hn 13 


-1.58 


0.05 


-1.04 


0.06 


-0.40 


0.05 


XX Cha 


-1.48 


0.04 


-0.81 


0.06 


-0.38 


0.07 


ISO 282 


-1.16 


0.08 


-1.48 


0.10 


-0.54 


0.07 
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Table 4 — Continued 



Object name 


"2-6 




"6-13 


C"n6-i3 


ril3-31 




T50 


-1.85 


0.10 


-0.98 


0.14 


-0.33 


0.11 


T51 


-1.09 


0.04 


-1.08 


0.07 


-1.52 


0.11 


CV Cha 


-1.41 


0.04 


-0.69 


0.04 


-0.27 


0.07 


T54 


-2.88 


0.07 


-2.38 


0.17 


1.10 


0.16 


Hn 21 W 


-1.98 


0.08 


-1.17 


0.09 


-0.77 


0.10 


T56 


-2.13 


0.12 


-1.14 


0.15 


0.90 


0.07 


2M J11241186-7630425 


-2.33 


0.08 


-0.81 


0.13 


0.68 


0.11 



2MASS K5 flux is an upper limit, therefore n2-6 value is not listed. 



Table 5. Properties of the Silicate features of Class II objects 



Object name Wio '^Wio ''Pio ^'ll.s/^'g.S "'Fll.s/Fg.s ^^^0 '^1*20 ^^0 "^^20 '^^3 "'W^^ 

(/im) (m™) (erg cm~^) (erg CTn~^) (m™) (m"^) (erg cm~^) (erg cm~^) (f^ra) 



SX Cha 


2.46 


0.08 


3.0B-11 


6.6E-13 


0.49 


0.03 


1.69 


0.23 


8.4E-12 


6.6E-13 


0.09 


0.06 


T5 


1.25 


0.07 


2.7E-12 


1.3E-13 


0.74 


0.07 


2.28 


0.31 


1.7E-12 


1.3E-13 


0.12 


0.08 


2M J10580597-7711501 


2.86 


0.15 


2.4E-13 


l.lE-14 


0.79 


0.05 


4.37 


0.37 


2.1E-13 


l.lE-14 


0.70 


0.11 


SZ Cha 


3.28 


0.09 


1.4B-11 


2.5E-13 


0.61 


0.03 


4.79 


0.35 


2.8E-11 


2.5E-13 


-0.01 


0.07 


TW Cha 


5.47 


0.12 


1.7E-11 


1.8E-13 


0.47 


0.02 


3.76 


0.34 


5.6E-12 


1.8E-13 


-0.08 


0.08 


OR Cha 


4.49 


0.11 


B.OB-11 


6.4E-13 


0.44 


0.02 


2.77 


0.2B 


1.8E-11 


6.4E-13 


0.04 


0.06 


CS Cha 


3.19 


0.09 


3.1E-12 


6.4E-14 


0.61 


0.03 


3.07 


0.32 


8.1E-12 


6.4E-14 


-0.01 


0.07 


CT Cha 


2.00 


0.08 


2.0E-11 


5.1E-13 


0.42 


0.03 


2.95 


0.26 


l.lE-11 


5.1E-13 


-0.00 


0.07 


ISO 52 


2.16 


0.08 


l.lE-12 


3.0E-14 


0.55 


0.04 


2.51 


0.32 


8.0E-13 


3.0E-14 


0.27 


0.08 



2M J 11062942-7724586 



CHSM 7869 


1.39 


0.08 


8.4E-14 


4.6E-15 


0.76 


0.06 


2.31 


0.38 


4.3E-14 


4.6E-15 


-0.14 


0.36 


ISO 79 


1.25 


0.07 


3.3E-13 


1.5E-14 


0.61 


0.06 


2.52 


0.32 


2.9E-13 


1.5E-14 


-0.13 


0.08 


Hn 5 


2.09 


0.08 


6.6E-12 


1.8E-13 


0.66 


0.04 


2.71 


0.32 


1.7E-12 


1.8E-13 


0.09 


0.10 


CHXR 20 


3.63 


0.10 


l.lE-11 


1.8E-13 


0.53 


0.03 


3.62 


0.34 


3.4E-12 


1.8E-13 


0.04 


0.08 


UY Cha 


2.95 


0.09 


4.0E-12 


7.8E-14 


0.65 


0.04 


2.17 


0.30 


9.0E-13 


7.8E-14 


0.40 


0.10 


2M ,711065939-7530559 


1.55 


0.07 


l.OE-13 


4.0E-15 


0.63 


0.06 


3.09 


0.33 


8.0E-14 


4.0E-15 


0.19 


0.08 



2M ,711070369-7724307 



ISO 91 6.93 0.14 5.1E-12 4.4E-14 0.35 0.01 4.76 0.35 3.1E-12 4.4E-14 0.03 0.07 

UZ Cha 2.15 0.08 4.3E-12 1.2E-13 0.61 0.04 3.82 0.36 2.6E-12 1.2E-13 -0.14 q.07 

CHXR 22E ■■■ •■ ■■■ ■■■ ■ ■■ ■■■ ■■■ •■■ ••■ ■ ■• ••■ 



Cha Hal 


2.01 


0.08 


5.3E-13 


1.6E-14 


0.58 


0.04 


1.95 


0.31 


3.2E-13 


1.6E-14 


-0.03 


m7 


T25 














3.39 


0.33 


1.6E-12 


2.2E-14 


0.15 


q.07 


DI Cha 


1.87 


0.08 


9.2E-11 


2.7E-12 


0.60 


0.04 


1.98 


0.24 


3.6E-11 


2.7E-12 


0.01 


0.06 


VV Cha 


2.29 


0.08 


5.0E-12 


1.2E-13 


0.45 


0.03 


2.05 


0.31 


1.6E-12 


1.2E-13 


0.06 


0.08 


Cha Ha2 


1.08 


0.25 


8.4E-13 


1.9E-13 


0.72 


0.31 


1.83 


O.Bl 


3.BE-13 


1.9E-13 


0.16 


0.26 


T28 


2.07 


0.08 


1.6E-11 


4.2E-13 


0.51 


0.04 


1.56 


0.29 


4.2B-12 


4.2E-13 


0.25 


0.08 


CHSM 10862 














1.68 


0.30 


l.lB-13 


l.lE-14 


0.14 


0.10 


CHXR SOB 


2.88 


0.09 


1.4E-11 


2.9E-13 


0.61 


0.03 


3.29 


0.32 


6.5E-12 


2.9B-13 


0.23 


0.08 


T29 


1.74 


0.07 


1.6E-10 


5.0E-12 


0.49 


0.04 


2.41 


0.24 


9.7E-11 


5.0E-12 


-0.01 


0.06 


CHXR 30A 


2.25 


0.08 


4.7E-12 


1.2E-13 


0.40 


0.03 


2.19 


0.30 


1.5E-12 


1.2B-13 


0.04 


0.07 


VW Cha 


2.42 


0.08 


6.6E-11 
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Table 5 — Continued 
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Note. — Values are not listed for objects for which a satisfactory continuum fit could not be obtained 
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Fig. 1. — Spectral type distribution of the Cha I sample for stars with spectral types later than 
KO. 








2 3 4 

Aj (mag) 



Fig. 2. — Distribution of extinction (Aj) towards the young stars in Cha I 
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Fig. 3. — SEDs of Class I objects in the Chamaeleon I sample. Photometric (DENIS, 2MASS, 
IRAC and MIPS) measurements are shown as solid cricles and the IRS spectrum as solid line. The 
dashed line represents the stellar photosphere. 
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Fig. 4. — SEDs of Flat spectrum sources in the Chamaeleon I sample. Symbols and line types have 
the same meaning as in Fig. [3l 
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Fig. 5. — SEDs of Class II objects in the Chamaeleon I sample, 
same meaning as in Fig. [3l 
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Fig. 5. — Continued. 
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Fig. 5. — Continued. 
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Fig. 5. — Continued. 
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Fig. 5. — Continued. 
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Fig. 6. — SEDs of transitional disks in the Chamaeleon I sample. Symbols and line types have the 
same meaning as in Fig. [H 
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Fig. 7. — SEDs of Class III objects in the Chamaeleon I sample. Symbols and line types have the 
same meaning as in Fig. [H 
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Fig. 8. — SEDs of stars which are probable members of e Cha group. Symbols and line types have 
the same meaning as in Fig. 21 
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Fig. 9. — Comparison of the observed spectral indices n2-25 and 115-12- Dashed horizontal lines 
indicate the regions occupied by Class I, Flat spectrum, Class II and Class III sources. Dotted 
vertical lines separate envelope, disk and photospheric sources. 
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Fig. 10. — IRS spectra of the Class I objects in our sample. For T14a, the dereddened spectrum is 
shown. Spectra have been scaled to fit the figure. 
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Fig. 11. — IRS spectra of the Flat spectrum sources in our sample. Objects for which Aj values 
are listed in Table 2, the dereddened spectra are shown. Spectra have been scaled to fit the figure. 
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Fig. 12. — Dereddened IRS spectra of Class II objects which show enhanced 10 fim sihcate emission 

(Wio > ^ 10, upper octile)- 
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Fig. 13. — Dereddened IRS spectra of Class II objects with moderate 10 fxm feature strength and 
flatter continuum 
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Fig. 14. — Dereddened IRS spectra of Class II objects with moderate 10 /xm feature strength and 
steeper (more negative) continuum slope 
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Fig. 15. — Dereddened IRS spectra of Class II objects with weak 10 fim silicate emission and flatter 
continuum 
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Fig. 16. — Dereddened IRS spectra of Class II objects with weak 10 fim silicate emission and 
steeper (more negative) continuum 
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Fig. 17. — Dereddened IRS spectra of Class II objects which could not fit in to the other groups. 
The first two are B type stars in the sample and the next two are outwardly truncated disks. 



-66- 




Fig. 18. — IRS spectra of transitional disk candidates 
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Fig. 19. — Continuum spectral index ni3„3i plotted against the index ng-is (see text for definitions) 
for Class II (open diamonds) and Class III (solid gray circles) objects. Transitional disk candidates 
are shown as solid diamonds. The star symbol represents location of the stellar photosphere {n = 
—3). The dashed line indicate ni3_3i = —4/3 appropriate for a geometrically thin, optically thick, 
fiat disk. 
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Fig. 20. — ni3_3i versus ng-is for Cha I objects. The symbols have the same meaning as in 
Figure [m Overplotted are the spectral indices computed from accretion disks models (see text for 
details) for 0.2 Mq (solid star symbol) and 0.5 Mq stars (solid circles). The model indices computed 
for accretion rates of 10^'^ MQyr^-*^ (top left panel), 10~® MQyr"-*^ (top right panel), 10~^ M0yr~^ 
(bottom left panel) and 10~^*^ M0yr~^ (bottom right panel) and for inclination angles of 60° (green), 
40° (blue) and 20° (red) are shown. Model indices for the settling parameter e of 1, 0.1, 0.01 and 
0.001 are represented by the sizes of the solid symbols from the smallest to the largest. 
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Fig. 21. — Distribution of continuum indices for the Class II (full) disks in Cha I (gray solid 
histogram) and Taurus (line histogram). 
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Fig. 22. — Continuum spectral indices n2_6 and ni3_3i plotted against each other for Class II 
objects in Taurus (top) and Cha I (bottom). The dotted lines represent the median of the indices 
and the dashed lines represent lower and higher octiles for the two regions. Transitional disks in 
both regions are shown as solid star symbols. 
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Fig. 23. — Equivalent width of the 10 silicate feature, Wio, plotted against the continuum 
spectral index ni3„3i for the Class II objects in Cha I. TDs and PTDs are shown as solid black star 
symbols and outwardly truncated disks are shown as solid black circles. The dotted lines represent 
the median in Wio and ni3_3i values and the dashed line the upper octile of Wio distribution. The 
solid gray circles represents objects with enhanced 10 silicate emission (Wio ^ ^ lo, upper octUe)- 
The polygon defines the region in the Wio - ni3_3i plane, allowed by the irradiated, accretion disk 
models for typical range of accretion rates, stellar masses, inclination angles and settling parameters 
e (see text). Typical errors are shown on the top-right corner (see Table [5]). 
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Fig. 24. — Integrated flux of the 10 fim silicate feature Fio plotted against the equivalent width 
Wio . The symbols have same meaning as in Figure [23l Typical errors are shown on the top-right 
corner (see Table Ej). 
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Fig. 25. — Equivalent width of 20 /xm silicate feature W20 plotted against the continuum index 
ni3_3i. The symbols have same meaning as in Figure [23l Typical errors are shown on the top-right 
corner (see Tabled]). 
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Fig. 26. — Ratio of the integrated flux of 20 fj,m silicate feature to that of the 10 /im silicate feature 
F20/F10 versus the continuum index ni3„3i. The symbols have same meaning as in Figure [23l 
Typical errors are shown on the top-left corner (see Tabled]). 
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Fig. 27. — Fu s/Fq^s plotted against the equivalent width of the 10 /im silicate feature Wio- The 
open star symbols represent TDs and PTDs in Tau-Aur region. The other symbols have same 
meaning as in Figure [23l SZ Cha is not shown here as it has a strong PAH feature at 11.3 fim and 
a reliable estimate of -F11.3/-F9.8 could not be obtained. Typical errors are shown on the top-right 
corner (see Table [5]). 
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Fig. 28. — -F'li.s/i^g.s plotted against the equivalent width of the 33 /nm silicate feature W33 for 
Class II disks in Cha I. The other symbols have the same meaning as in Figure [23l SZ Cha is not 
shown here as it has a strong PAH feature at 11.3 /im and a reliable estimate of -F11.3/-F9.8 could 
not be obtained. Typical errors are shown on the top- left corner (see Table E]). 
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Fig. 29. — -F'li.s/i^g.s plotted against W33 for Class II disks in Taurus. The solid star symbols 
represent TDs and PTDs. Typical errors are shown on the top-left corner (see Table [5]). 
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Fig. 30. — Median SEDs of Class II objects in the spectral type range of K5-M2 (top panel) and 
M3-M8 (bottom panel). The gray shaded regions are bound by the lower and upper quartiles of the 
median in both the panels. Synthetic SEDs generated from accretion disk models for a 0.5 Mq star 
(top panel) for accretion rates of 10^''^ (blue), 10~^ (red) and 10"^ (green) MQyr~^ and for depletion 
factors e = 0.001 (solid line), e = 0.01 (dashed line) and e = 0.1 (dotted line) are overplotted on 
the K5-M2 median. SEDs generated from accretion disk models for a 0.2 Mq star (bottom panel) 
for accretion rates of 10~® (red), 10~^ (green) and lO"^'^ (cyan) M0yr~^ and for depletion factors 
e = 0.001 (solid line), e = 0.01 (dashed line) and e = 0.1 (dotted line) are shown for the M3-M8 
median. The models shown are for an inclination angle of i = 60°. The photospheres plotted (solid 
black line) are for a K 7 star (K5-M2; top panel) a nd for a M5 star (M3-M8; bottom panel). The 
photospheres are from lKenyon HartmannI (|l995l ) and are normalized at J-band. 



-79- 




Fig. 31. — Distribution of the continuum indices for close (< 1") and wide (> 1") binaries and 
single stars among the Class II objects in Cha I 
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Fig. 32. — Frequency distribution of the spectral indices n2-6 (a), ^^6-13 (b) and ni3_3i (c) for the CTTS (top panel) and WTTS 
(bottom panel) in our sample. The solid histogram represents TDs and PTDs among CTTS and WTTS. 



